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1. INTRODUCTION 

The aerosol backscatter coefficient 0 * m^sr' 1 ] 

is defined as the fraction of incident electromagnetic energy at 
wavelength A [m] which is scattered by atmospheric aerosol 
particles in the backward direction per unit solid angle, per 
unit atmospheric length (Measures, 1984, pgs. 53 and 239, and 

Kavaya and Menzies, 1985). It is assumed to be constant over the 

measurement volume of interest. In general, the atmospheric 
scattering (including backscatter ing) will be due to a complex 
distribution of scatterers that vary in composition, size, and 
shape. In order to account for arbitrary polarizations of the 
incident and scattered (detected) light, 0 must be generalized to 
a 4 x 4 Mueller matrix (Kavaya, 1987; Anderson, 1989) . Once the 
Stokes vectors of the incident and detected light are specified, 
e.g. for a specific laser radar (lidar) system, then the 
backscattered energy is proportional to a weighted sum of a 
subset of the 16 0 Mueller matrix elements. The weighted sum may 

be considered an effective scalar value of 0. This will be 


assumed here. 

A Light Detection And Ranging (lidar) system may measure 0 
in many different ways: 1) the lidar may be continuous -wave (CW) 
or pulsed; 2) the lidar may employ direct detection or heterodyne 
detection; 3) the laser beam may be focused or collimated; and 4) 
the measurement volume may contain numerous aerosol particles 
(volume mode - VM) , a few particles, or only one particle at a 
time (single particle mode - SPM) . The NASA GLObal Backscatter 
Experiment (GLOBE) program is attempting to increase the 
knowledge data base of global aerosols and their backscatter 
coefficient with particular emphasis on assisting the design of 
the proposed NASA Laser Atmospheric Wind Sounder (LAWS) lidar 
system. This Earth-orbiting wind profiling lidar will depend on 
aerosol backscatter for its signal. It will most likely be a 
pulsed, heterodyne (coherent) detection, collimated lidar system 
making global 0 measurements in the volume mode over the entire 
troposphere. If it orbits at a height of 800 km and if the laser 
beam has a nadir angle of 52.7°, for example, then the 
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interaction of the laser light with tropospheric aerosol 
particles will occur at a range of 1500 km with an effective 
(local) nadir angle of 63.5°. The relative Doppler shift between 
the laser light and the aerosol particles will depend on the 
satellite velocity, the lidar's azimuth scan angle, the latitude, 
and the wind velocity. 

This contract effort addressed the measurement and 
calibration of 0 by a CW, coherent detection, focused lidar 
system operating either in single particle mode or in volume 
mode. It will be shown that within reasonable assumptions, the 
value of 0 measured by the CW VM or SPM lidar is identical to the 
0 value which affects the LAWS instrument performance. 

Other important issues of position (latitude, longitude, 
a lti.tude) t spatial resolution, time, temporal resolution, 
Doppler shift, incident and detected polarizations, 
atmospheric conditions, origin of air mass, recent weather 
conditions, measurement volume size and shape, and lidar 
beam-to-measurement volume orientation are not addressed here. 

In Section 9 we describe the work performed in development 
and testing of a 2-/ra Tm , Ho : YAG coherent lidar system. The 
system was recently moved to the Table Mountain remote field site 
(~ 10 miles north of Boulder, CO) . Preliminary hard target and 
atmospheric data are presented in Section 9.3. We believe this 
to be the first demonstration of a coherent eyesafe solid-state 
lidar system. 
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2. OVERVIEW OF THEORY 


The aerosol backscatter coefficient P is a measure of the average power backscattered by the 
atmosphere. For coherent detection laser radar, the transmitted beam is confined to a narrow angular 
region of the atmosphere. This permits the propagation, scattering, and reception of the fields to be 
described by Fresnel diffraction theory. For coherent detection lidar, the intermediate frequency 
(IF) signal from the detector I s (t)[A ] is the basic signal for generating statistics to estimate p. Since 
the IF signal is a random quantity, a statistical estimate for P (a function of IF signal samples that 
approximates the true value of P ) is required. 

A satisfactory estimate for P requires that it is an unbiased estimate of P and that the estimate is 
accurate. Both volume mode (VM) and single particle mode (SPM) estimates of P are unbiased for 
ideal operating conditions. The ideal operating condition for the SPM method requires that the 
"threshold level" for identifying single particle events is low enough to sample enough of the 


scattering particles that contribute to p. 

The accuracy of the estimates is affected by three main sources of error: 1) the additive noise 

of the detector, 2) the random fluctuations of the heterodyne power, and 3) the systematic errors of 

. . ^ L AP (noise) AP (random) 

the system. The relative errors of these three components are denoted by ^ , P ’ 

and ^systematic)-' respect ively. If all these errors are small, then the total relative error (i.e., 

P 

accuracy of estimate or figure of merit) of the estimate is 


P 


[ A$(noise ) 

2 

f * 

AP( random ) 

2 

r ' 

^(systematic ) 

J 

P 

L v J 

T 

P , 


0 J 



( 1 ) 


Note that AP (noise) is related to the signal-to-noisc ratio (SNR), AP(random) arises from the ran- 
dom heterodyne detection process, and AP(systematic) is linked to sensitivity analysis. 

The merit of any unbiased estimate of P is then given by this total relative error. The first two 
error terms will be considered in this report, for both VM and SPM operation. The systematic error 
requires involved calculations of the actual lidar system parameters. This analysis may be the most 
important since systematic errors may be the dominant source of error. Sensitivity analysis is one 
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part of systematic error analysis. Quantitative analysis of performance will be obtained using the 
current Marshall Space Flight Center (MSFC) CW coherent lidar system parameters [W.D. Jones, 
personal communication, 12/9/89 & 1/4/90]. 


3. VOLUME MODE (VM) MEASUREMENTS OF 0 


3.1 Introduction 

A coherent detection lidar [See Figure 1] consists of a transmitter laser described by the scalar 
field Ej (tf,z )[(Wm 2 ) 1/2 ] where If = (u. My )[m ] is the transverse coordinate at a propagation distance 
z[m]; a dimensionless effective receiver lens JV(V) where V= (v x ,v y )[m] is the transverse coordi- 
nate in the receiver plane z — 0; a detector in the plane z — L described by the quantum efficiency 
function T\(&)[electrons /photon] at transverse coordinate W[/n] on the detector surface; and a 
local-oscillator (LO) field li£/?(V,z). The irradiance profile of the transmitter beam Ij(\i,z)\Wm~^} 
is given by 

h( fc) = E t ($j)E t ’ (tf^z) (2) 

and normalized such that 


P T - (3) 

where P T [W] is the power of the transmitted beam in the absence of extinction. 

An analytic expression for lidar performance can be obtained for a Gaussian lidar system; i.e., 
when 


is the transmitter field, 


£Y<ttO) = 


P T 

1/2 

2 

exp 

Krf 





(4) 


4 


W(V) = exp 


(5) 


v 2 



ikv 2 
2 F r 


is the receiver lens, and 


Euo CM) = 


P LO 

1/2 

w 

v 2 

iky 2. 

X 
1 

exp 

hi 

2 rjj) 

‘ 2F lo \ 


( 6 ) 


is the LO field backpropagated to the receiver plane z =0. The parameters r T ,r R , and r w [m] 
represent e~ l (37%) intensity radii; and the parameters Fj,F r , and F^q [m] represent phase curva- 
tures. (The actual MSFC CW coherent focused lidar systems may have truncated Gaussian transmit- 
ted and LO beams due to a 15 cm diameter rectangular or "top hat" receiver mirror function. The 
solution for this case is much harder and may not lead to an analytic expression.) The current MSFC 
CW lidar system will be modeled as a matched monostatic Gaussian lidar 


r T ~ r LO ~ 0-02157 m , transmitter and local oscillator 1/e intensity radii (1/e 2 radii = 0.0305m) 

r R = oo, no truncation of the LO by the receiver 

F T =F r =F = 10.59 m , transmitter and receiver radius of curvature 

F lo = oo, local oscillator phase radius of curvature at receiver plane 

\ = 9.114519 pm ,C0 2 laser wavelength [9R(20) line of 12 C 18 0 2 isotope] 

= 0.12, optical transmission + detector AC quantum efficiency [electrons/photon] 

P T = 4 w , transmitted power 

B vm = 100 kHz , VM bandwidth (SAW) 

B SPU = 1.5 MHz , SPM bandwidth 
V i4 7cos0= 100 m 
V A [m/s] = airplane velocity 
T[s] = Is observation time 

0 [rad] = angle between lidar optical axis and the perpendicular to the airplane velocity vector 
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Volume mode operation occurs when many scattering aerosols are illuminated by the 
transmitter laser. The statistics of the backscattered field is then described as a complex Gaussian 

I 

random process which is sometimes called speckle statistics. The backscattered field from the aero- 
sols is collected by the receiver lens and mixed with the local oscillator field in the detector plane. 

The square of the IF signal ig [A 2 ] (obtained using a square-law-detector) is the coherent IF power. 

If the detector is larger than the local oscillator beam and if the detector quantum efficiency is uni- 
form over the illuminated area of the local oscillator, and if the additive shot noise is negligible, then 
the average IF power is related to P by (see Appendix A) 

<//> = p G v P t (7) 

where <•> denotes ensemble average, G v is the system gain [A 2 m sr IT -1 ], i.e., 

G v = $G v (R)dR (8) 

o 

and 

* >2 mm mm 90 

Cv(S) ' 2 'T v I I I <W)co' l #)cV 1 ,W)> 

W^)W* (^ 2 )E lo * (*„o: )E lo (V 2 ,0)dV,dV 2 d^ (9) 

is the volume mode coherent lidar system range gain density [A 2 sr W -1 ] as a function of range 
/?[m]. Here, G D = dimensionless amplifier gain, h = 6.626x1 0 -34 [/s] is Plancks constant, and 
GOWOtm' 2 ] is the Green’s function; i.e., the field at (V,0) due to a unit point source at (j fjt). 
Under the Fresnel approximation and no extinction 

G(W) “**®‘ P (10) 

where k = 2ji fk [m _1 ] is the wavenumber of the optical field. The system range gain density func- 
tion identifies the dominant scattering regions as a function of range. 
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The aerosol backscatter coefficient |3[m 1 sr -1 ] is given by 

oo 

P = Jon(o)do (11) 

0 

where n(a) [m'V] is the number of aerosol particles per unit volume per unit cross-section 
o[m 2 jr _1 ] and is assumed to be constant over the measurement volume defined by G V (R). The 
system gain G v can be determined from the system components which requires integration of Eq. 
(9), which is difficult for a real system. Alternatively, the system gain G v can be characterized by a 
measurement of Gy(R ) at appropriate ranges R. 

Since a measurement of lidar IF power includes the additive shot noise signal i n (r), the correct 
statistic for VM measurement of P is 

P -P (signal )-P ( noise ) (12) 

where 

i r 

P (signal ) = ^rfi 5 2 (r )dt (13) 

1 o 

is an estimate of the IF coherent signal power [A 2 ] with shot noise included for an observation time 
T[j], and 

i T 

P (noise) = MQ(t)dt (14) 

1 0 

is an estimate of the IF noise power [A \ The estimate P [A 2 ] is an unbiased estimate of P, i.e., 

<P > = P G v P T . (15) 

The dimensionless lidar SNR is the average coherent lidar power <Q> normalized by the aver- 
age noise power <t„ 2 >. For the given assumptions 

SNR =$G snr P t (16) 

where the system gain [m sr W 1 ] for SNR is 
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°SNR = T vW ~J J / K(WX?^JW)C*(W)> 

nvorix) o - - . » . - 

w OfyW'* * <?1 S»Euo ^O^id^^dR (17) 

where 

Plo = ] \Eu,(te)\ 2 dt (18) 

is the LO power [W] measured by the detector. Then 

g snr - fGsNR(P)dR > (19) 

o 

where Gsnr (R ) t* r W~ l ] is the SNR range density function. G$nr (R ) can also be determined from 
system parameters but this is difficult to do accurately. 

The system range density function (Gy (R ),G$nr (R )) indicate the regions of propagation range I 
that are most important for system performance. Another measure of this behavior is the cumulative 
distribution function (CDF) defined by 

R 

\g v (r ym' 

CDF (R ) = (20) 

\G v (R)dR 

o 

which indicates the fraction of the signal statistic from aerosols with range less than R . 

For the Gaussian lidar system, the range density function is 



where 
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( 22 ) 


_L = _L + _L- 

ri r R r LO 


and K(R ) is the dimensionless one-way irradiance extinction. 

For the Gaussian lidar system, the SNR range density function is 

4iO\r£[K(R)] 2 


1 SNR 


(*) = 


hvBvM r LbR 2 


11.,. R Z- + S-? 


(23) 


In the limit as R ->0, the range density function approaches a constant Normalizing the range den- 
sity function by this constant produces 

k\rf + rg) 


Gytf(R ) - 


1 . 1 


. * R , n R I S. 

— X" ■*" o 0 t? ) - o ^ V I7_ _ d 2 


2-2 


rf ' r £ 


/?' 


^ Flo R- 


(24) 


Performing the integrals in Eqs. (8) and (19) for the matched monostatic model of the MSFC lidar, 
uniform P, and constant extinction K produces 


Gy — 2XP jjQ 


KG D r\e 

2 

£ + tzn\kr 2 /F) 

h\ 


2 


(25) 


and 


Gsnr - 


ML iL + tan ~\kr 2 /F) 
hvB VM 2 


(26) 


The total system gain for a focused lidar will vary by only a factor of two for any change in radius of 
curvature F or beam dimensions r , assuming uniform P and constant K. Therefore, these parame- 
ters can be chosen to improve other aspects of system performance. 
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3.2 Error Due to Additive Shot Noise 


For the VM estimate of P, a measure of the shot noise component of error is defined as 

Ap (noise ) {VAR [P ( noise )1) 1/2 <Q> 

P <P> <P>{TB Vht ) m 


( 27 ) 


where B VM [Hz ] is the VM bandwidth of the square-law-detector and the noise power has been accu- 
rately estimated from a calibration run. For the MSFC matched monostatic system using the SAW 


spectrum analyzer with a 100 kHz bandwidth, 


AP (noise) 5.02xl0~ 13 

0 n\$P T T\cosQ(TB VM ) m PcosO 


(28) 


For cos0 = 1 and P = 10~ 12 m~ l sr~ l , . 4Pfo > i,se ) _ o. 5 t which is detectable (50% error). An observa- 
tion time of T = 10 s would produce - =0.16, and an observation time of T — 100 s 

would yield - ^ n ^ ise I = 0.05 (5% error). 

P 


3.3 Error Due to the Random Process 

The error of the VM estimate for P due to the random fluctuations in heterodyne power is the 
dominant source of error for high SNR and negligible systematic error. For ideal VM performance, 
the heterodyne power has "speckle statistics," i.e., the IF current is a zero-mean complex Gaussian 
random process and the heterodyne power has a negative exponential distribution. For focused VM 
operation, the IF power may be characterized by single particle events. The IF power then has non- 
ideal statistics (not speckle statistics) and a more complex analysis must be performed. 

The variance of the VM estimate P [Eq. (12)] due to speckle statistics is 

<P > 2 T d 

VAR IP] y-i- (29) 

where T d [r ] is the temporal decorrelation time of the "speckle process." The accuracy of P is then 
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A Qfaindotn ) _ (VA/? [P ]) rj» ^p-ii/2 

P " <P> d 


For T d = liu, N^randomX _ q. 001 = 0.1%, assuming a T = 1 s observation time. 

P 


( 30 ) 


3.4 Calibration of VM Measurement of P 

VM measurements of P can be calibrated using a diffuse target with reflection coefficient 
p [sr -1 ] . The backscattered irradiance I (?, 0) [Wm~ 2 ] at the receiver is given by 

nrn-^f < 31 > 

R 2 

where Pq [W] is the calibration power transmitted onto the target. The average coherent IF power 
from the hard target at range R is then 

<//(/? )> c =p/ > c G v (/?) (32) 


yielding 


G V (R) = 


<i S \R ) > c 
P p c 


as a measure of G v (R ). The system gain G v is obtained by integration of G V (R ) with 
As with calibration of lidar power range density, a measurement of SNR from 
diffuse target is a measure of Gsnr ); i e * 


(33) 

respect to R . 
a known hard 


Gsnr ( r ) - 


SNR (R ) 

P p c 


(34) 
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4. SINGLE PARTICLE MODE (SPM) MEASUREMENTS OF p 


4.1 Introduction 

A focused CW lidar beam will have an effective sensitivity or gain along its optical axis which 
is peaked at the focus. For a tightly focused lidar beam, and for sufficiently low values of (3 (this 
implies low particle concentration), the received signal power will be characterized by randomly 
occurring events as particles of different composition, size, and shape traverse the focused beam at 
different ranges along the optical axis R [m] and at different distances p x [m ] above or below the 
optical axis. We assume the lidar system is moving through the atmosphere with velocity V A in the 
direction of the + y axis (e.g., in an aircraft) and that it is aimed in the y-z plane at an angle to the 
+ z axis (see Figure 2). For example, 0 was + 17.44® during a recent flight measurement program 
(Gras et al, 1988). 

The m ax im u m lidar signal Sj(p x ,p ym ,R) [A 2 ] for the jth particle event will occur at some 
coordinate (p x ,p ym ,R) = (] t m ,R) of the trajectory, where p x is the displacement from the laser 
beam axis. The proper statistic for estimating (J in the absence of additive noise (see Section 4.6 for 
correction for noise spikes) is the sum S [A 2 ] of the maximum (or peak) signals: 

N 

Sj(P x .P ym >R) (35) 

7 = 1 

If the threshold for the identification of peak signals is low enough to identify all the scattering parti- 
cles, then 

<S > = TP t G s cosG (36) 

where G$ [A 2 W~ l sr] is the total system gain given by 

G s =]] G s (p xr R)dp x dR (37) 

0 - oo 

and 
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lilsiLii 


G s (p x Jl) = 2P T 


r ]G D e 


hv 


2 7 7 j T ® m X) 

1 j — — cx p 


«■ 


^[tt»-*i> 2 -W.-*2) 2 ] 


W&iW* tfdEuo • ^,0)£ w OkOVfVj^ 


( 38 ) 


is the system’s transverse area gain density [4 2 W~ l m 2 sr ] for a scattering particle at (jf OT *R ) and 


/ r (ft*) = 


Ij^R) 

P T 


(39) 


is the normalized transmitter irradiance [m~ 2 ] at transverse coordinate Jf [m ] and range R . For the 
general Gaussian lidar system, the ensemble average of the normalized transmitted beam profile is 
then 


<J T (pJO> = 


r > 



(40) 


where 


R 


rim - r T \ l-^) 2 + (^-) 2 + 2 [ W) 




(41) 


where r B is the e~ 1 intensity radius [m] of the ensemble averaged transmitted beam profile, and 

1 - 3/5 

I A „ p' 

Sq(R) s 


Hk 2 jC 2 (R')(l--^)dR' 
o K 


(42) 


is the field coherence length [m] due to refractive turbulence, H = 2.914383, and C 2 (R ) [m -2/3 ] is 
the refractive index structure constant profile. The last term in Eq. (41) is the spread of the transmit- 
ted beam profile due to Kolmogorov refractive turbulence with zero inner scale. This effect must be 
negligible for simple calibration. If it is not negligible, then the calibration will depend on refractive 
turbulence conditions. We assume it is negligible in this paper. 

The transverse area gain density function G s (p x Jt) for this idealized system is 
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Gs(p x J*) = 8 


ri G D e 

P W r E 

2 2 
Px Px 

h\ 

<. * 

k 2 rjj) r R (R )ry(R ) ° XP 

rj(R) r$(R) 

* * 


where 


-L = _L + JL_ 

r i r R r L0 


and the effective receiver area is Wg , and 


'«(* > = ri (1 - * -2-f + (-f-) 2 

t R r LO kr E 


(43) 


(44) 


(45) 


r M (R) [m] is the contribution to the gain function from phase mismatch between the received field 
and the local oscillator. 


The range gain density function G s (R )[A 2 W~ 1 m~ l sr ] is obtained by integrating G s (p x Jt) 
overp x ,i.e. 

^2 - .2 


G s (R) = $k 


1/2 


r\G D e 


hv 


P LO r E 


k 2 r& \ri(R )r$(R ){rj(R ) + rJ(R ))] 


1/2 


(46) 


It is highly recommended that both theoretical and experimental sensitivity analyses of G s (R ) 
be performed about a candidate operating point of the lidar system. 

In the limit of small R 


G s ( 0) = 8it 


1/2 


TlGpg 

hv 


P LO r E 


k lr wl r T r i( r T + r £)] 


2m1/2 


(47) 


G S (R ) becomes independent of R for small R just as the VM range gain density function G V (R ), 
did. We may normalize Eq. (46) with Eq. (47) to obtain 


Gsijt ) [ r T r E ( r T + r £)] V2 

G S (0) ~ [r£(R )r& (R )(r B 2 (R )+r£(R ))] 1/2 


(48) 


For the MSFC matched monostatic system, [r^ =r T = r , F w -r R * <*», F T = F R =F , 

c„ 2 =0] 
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(49) 


./»>-»* 4* + (£>* 


and 


2 2 2 
r E = r L0 =r 


(50) 


Assuming/? =F 


k 2 r 6 

GaW-JS^i 


(51) 


If/? 


Gsh(R 


r 3 (3/4) 


1/2 


W 


For tightly focused conditions, it is safe to assume F « Jr 2 / 2, so that Eq. (52) becomes 

f 3 

Gcw (/? — >°°) = — 7 


/?■ 


(52) 


(53) 


If F » kr 2 , Eq. (5 1) should be used with R substituted for F . 


4.2 Error Due to Additive Shot Noise 

For an unbiased estimate of P in SPM operation, the threshold level for signal spike 
identification must be low enough to sample the majority of aerosol particles that contribute to (J. 
When this is satisfied, the contribution of noise spikes that are erroneously counted as signal spikes 
is negligible unless (3 is low and few spikes are observed. The error analysis will assume that the 
number of signal spikes is much larger than the number of noise spikes. 

For the SPM estimate, the error due to additive shot noise is 

A$(noise ) _ (VAR [S ( noise )]) 1/2 _ > > (54) 

P <S> " <S> <N spike > m <S j > 

where <N spike > is the average number of spikes observed. For the MSFC matched monostatic Gaus- 
sian system 
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The function H(p x Jl) is the normalized system gain at range R and offset p x . The threshold is 
defined as e</ n 2 > and it is assumed that the contribution of the noise spikes above the threshold is 
negligible (see Section 4.6 for noise spike correction). For useful operation, e>3, provided there are 
more signal spikes than noise spikes. If e is not high enough, especially for low (3, the number of 
noise spikes will become larger than the number of signal spikes, and it will be difficult to accu- 
rately correct for the noise spikes. For a matched monostatic Gaussian system with no beam trunca- 
tion 

TTJ' X P - 2 7T 

k r r B 

The average number of peak signals above the threshold level depends on system geometry and the 
form of the distribution of scattering cross section n (a). 
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43 Error Due to the Random Process 

The accuracy of the SPM estimate of 0, due to the random fluctuations of heterodyne power, is 
determined by the fluctuations in peak signal levels, i.e.. 


A ^(random ) _ (VAR[S]) 1/2 _ 
(J <5 > 



•o *e 

1/2 

h m 

J \H 2 (p x Ji)dRdp x 
--0 



(60) 


(3[V A T cos0] 1/2 J J H(p x Jt)dRdp x 
-«o 


where 


^2 = j<7 Z n(c)do 
0 


(61) 


1/2 


^2 

is the second moment [m rr -2 ] of particle cross section. A small ratio of — jp is desired (a narrow 
distribution of a values for a small SPM random error). For the MSFC matched monostatic Gaus- 
sian system 
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P[15V a Tcos0] 


An accurate estimate for P requires sampling many spikes, which samples the many different sized 
aerosols. If p is significantly affected by large particles, then the estimate for p can be poor since it 
would require a long time to observe the rare events due to the relatively few large particles. The 
first and second moments of P provide the insight into the importance of different sized particles 
within n (a) for accurately estimating P . 


4.4 Calibration of SPM Measurement of P 

The lidar system can be calibrated by shooting test particles with cross section c c through the 
laser beam of power Pq , at an angle - 0£ , beam offset p K , and range Rq . The particle s cross 
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section ct c may be determined theoretically or experimentally. The maximum (peak) lidar signal 

Sc 2 ] is 

s c (Pxc’ R c) = F c G c G s (P x J*c ) (63) 

Therefore, 

G s (Px -Re ) = Sq iPxc R c V (F c ) (64) 

may be experimentally determined. Ideally the measurement is deterministic, but averaging of mul- 
tiple test particle shots may be employed to reduce noise, trajectory jitter, etc. 

The total system gain G s requires a 2-D integration over range R c and offset p K . The integra- 
tion overp^ produces the system’s range gain density [A 2 W~ l m - l sr] as a function of range , i.e., 

O s (Rc) = I G s (Px R c )dPxc , (65) 

which allows identification of the region (in R) of the maximum contribution to the overall system 
gain G s . This region will be near the focus of the lidar system. 

The dependence of the functions G s (p xr R) and G s (R ) on various system parameters is impor- 
tant for optimal design of a CW SPM lidar system, and for optimum calibration methodology. The 
calibration apparatus should span values of p x and R sufficient to include, for example, all values of 
Gs (p x R ) down to 1 % of its maximum value. 

Assume that B c = 6 , i.e., the trajectory of the calibration test particles through the laser beam 
is at the same angle as the trajectories of atmospheric aerosol particles during flight. Then, an esti- 
mate for p is [see Appendix B, Eq. (B6)] 

N 

o _ j m 1 

p t g s v a t cose (66) 

The cos0 dependence reflects that the effective sensing volume of the lidar decreases with increasing 
® • For the same SPM statistic S , a larger value of 0 implies a larger value of p. 
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Table 1 s ummarizes the assumptions made in arriving at Eq. (66) for calibrated SPM. Appen- 
dix C discusses the calibration methodology if the particular value of for a test particle is not 
known, but the probability density function (PDF) of the values of p x for many particles is known 

to be uniform. 
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Table 1. Assumptions of Calibration Methodology 


1. Fresnel approximation (narrow angle scattering). 

2. Negligible refractive turbulence. 

3. Sufficient resolution and coverage in p K and R c during calibration to accurately determine sys- 
tem gain. 

4. 6 C = 0. 

5. The test particle’s velocity equals the atmospheric measurement lidar system velocity or the 
lidar system receiver gives the same peak signal reading for different particle speeds (i.e., flat 
bandwidth). 

6. Known value of a c of test particles. 

7. Good linear peak power detector. 

8. Lidar system unchanged between calibration and aerosol backscatter measurements. 

9. P represents weighted sum of 16 Mueller matrix elements. 

10. Statistical process yielding peak signals is ergodic and stationary. 

1 1. n (a) constant over dominant region of measurement vol ume , 

12. P uniform over measurement volume. 

13. Sufficiently large number of particle events during SPM atmospheric measurement. 

14. Sufficient SNR. 

15. The mean heterodyne signal power is available. 

16. The heterodyne detector is bigger than the LO beam. 
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43. Requirements for SPM Operation 
The requirements for accurate SPM estimates include: 

1) Heterodyne power signal with well defined peak signals. 

2) Enough peak signals above threshold to produce unbiased estimates of (3 . 

3) Errors due to shot noise and statistical fluctuations in peak signals require knowledge of n (o). 

4) Errors due to misali gnme nt and optical mismatch (beam size and curvature) require numerical 
calculations. 

5) Accurate method of obtaining peak signals: 

a) A© must be large enough to sample spike events. For the MSFC system and V A = 100 m/s , 
the spike events are 7.1pj (X*, in Figure 3). 

b) Either average IF power (dashed line in Figure 3) is used or peak IF power (solid line in Fig- 
ure 3) with no filtering is used. If x d «x w , then using a smoothing filter with time constant 
Xf , where x d < Xf <x w , will produce a more accurate signal where the peak value is the average 
IF power of the signal. If x d =t w , then the maximum IF power must be used, and then the 
peak power Pp ^ = 2<Ppeak >* where <P p€a t > is the peak of the average IF power <j$ 2 > 
shown as the dashed line. If x d > 1 * , the IF frequency is not high enough to sample the single 
spike event and SPM is not recommended. 


4.6. Noise Signal Correction for SPM 

For low p, the number of signal spikes from aerosol scatterers will be small and the contribu- 
tion from random noise spikes must be considered. An estimate that removes the contribution of 
noise spikes is 
N 

S = 2 ) SjiPxiPynJt) “ < ^spik£(noise)' > <^noise ^ spike ^ spike (noise ^ ($7) 

j * i 

where <N spike (noise y* is the average number of noise spikes above the threshold and <5^^ > is the 
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average lidar power of the noise spikes. The second term connects the contribution from noise spikes 
and the last term corrects the contribution of the additive noise power to the signal spikes. When the 
number of noise spikes is comparable to the number of signal spikes, the error analysis must be 
extended to include the noise spike contribution. 


4.7. VM Random Error in SPM Regime 

The VM estimate for (S is unbiased even when operating in the SPM regime (i.e., the measure- 
ment volume does not include enough particles to cause speckle statistics). However, the random 
fluctuations of the heterodyne power are not characterized by exponential statistics but determined 
by the random fluctuations of the single particle events. The accuracy of the VM estimate due to 
these random fluctuations is given by 


P 


P 2 l/2 

J jF 2 (p x R)dRdp x 

1/2 

(VAR [F]) 1/2 

-~0 



<P> 


Pt^Tcos e] V2 j jF(p x Jt)dRdp x 
-~o 


( 68 ) 


where 


7 7 7 'rak*) 

F <Px*) m J J J ex P 




---- - R^w 

W^)W* o *i>E w * <? x Wto (69) 

is the integrated IF power [W] for a single particle event with offset p x . For the tightly focused 
MSFC matched monostatic Gaussian lidar system (F « kr 2 ) 


A Random ) ^fo 172 

P " 7t 5/4 P[6FrV 4 T cos0] 1/2 


(70) 
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5. SYSTEMATIC ERRORS 

There are many systematic errors in lidar calibrations. These include nonuniformities in quan- 
tum efficiency over the detector; non-linearities in the photo detector and square-law-detector, errors 
in the parameters of the transmitter field, LO field, and receiver response; and fluctuations in all 
parameter settings after calibration is performed. The magnitude of the systematic errors that 
involve system par ame ters can be determined by a sensitivity analysis, which identifies which 
parameter settings cause the most change in system gain due to small changes from their nominal 
settings. A major source of sensitivity for coherent detection lidar is beam alignment Ideal opera- 
tion of coherent detection lidar requires an excellent match between the transmitted Gaussian beam 
and the back-propagated local oscillator (BPLO) beam, especially for tightly focused conditions. For 
VM operation and equal transmitter and BPLO radii, the effects of beam axis offset <f [m ] and point- 
ing angle error A(? [rad] are given by 

G V (R) = Gyj (R )exp[-(A $R + ~f$) 2 l(2rjf(R ))] (71) 

where GyjiR ) is the ideal system gain density with no misalignment between the transmitted beam 
and local oscillator beam. The effects of misalignment can be easily calculated. 

For the focused system, the major contribution to (3 comes from a small region around the 
focus. An angle error of 67 \i rad between a perfectly matched transmit beam and local oscillator 
beam would produce a 40% error in the estimate of P . An offset of 0.71 mm between a parallel 
transmit and LO beam would also produce a 40% error. 

The error in the SPM estimate of p due to misalignment of the transmitter and LO beams can 
be calculated in a similar fashion as for the VM case. The behavior will be more complex but have 
a similar scaling. 

SPM measurements suffer from a bias since a threshold is required to determine if a single par- 
ticle event occurred. A measure of this bias is the ratio of the ensemble average of the SPM statistic 
with a threshold <S rt > to the ensemble average with no threshold, i.e. 
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( 72 ) 


- „ 1 1 D„lo c (p,Jl)lH(p : ,Jl)<ip x dR 
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D N [x] = ± 

, P 

J an {a)d a 
o 

is the fraction of P due to aerosol particles with a cross section a larger than x . 
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6. COMPARISON OF VM AND SPM 

6.1 Comparison of Range Gain Density Functions 

Figure 4 shows the normalized range gain density functions for both SPM and VM plotted vs 
RIF for the MSFC matched monostatic Gaussian system. Both curves become independent of R 
and equal to 1 for small R . The SPM G s (R ) falls as R~ 3 for large R [see Eq. (53)] and the VM SNR 
falls as R~ 2 for large R [see Eq. (24)]. This difference is due to the effective 1-D integral through the 
Gaussian beam in SPM and the effective 2-D integral over the entire Gaussian beam in VM. For 
large R the beam waist is proportional to R . Therefore die 1-D integral is proportional to R~ l and 
the 2-D integral is independent of R . Both SPM and VM have a further R~ 2 dependence due to the 
decreasing solid angle of the receiver aperture. It is difficult to tell from Figure 4 how big an interval 
in range must be included to account for 90% of the gain, for example. This requires the cumulative 
density function (CDF) of G S (R) and G V (R) [see Eq. (20)] which is plotted in Figure 5. Since 
Ft = Fr = 10.59 m, the abscissa values of 0.9 and 1.1 correspond to R = 9.53 and 1 1.65 m, respec- 
tively. A slightly larger range interval is required for the VM case for equal percentages of the sig- 
nal. Figures 6 and 7 show the same plots with the larger focal range and smaller beam radius used 
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in a recent NASA GLOBE flight (W. D. Jones, personal communication, 4/20/90 and 9/18/90). 
Since F T = F R * 50 m in Fig. 7, the abscissa values of 0.5 and 1.5 correspond to R=25 and 75 m , 
respectively. As expected from Eq. (51), the SPM peak at R =F (=90) is approximately a factor of 
9 l(T 3 smaller than the SPM peak at R * F in Fig. 4 (=10 4 ). Figure 8 shows the SPM G$ (R ) and 
the VM Gy {R ) for the same conditions as Figure 4 except X = 1.064158 pm. The focusing is much 
tighter with the shorter wavelength. As expected from Eq. (51), the SPM peak at R = F (=6 10 6 ) is 
approximately a factor of 628 higher than the SPM peak in Fig. 4. The tighter focusing is also exhi- 
bited in the CDF curves of Figure 7. 


62 Error Comparison 


Calculations of SPM performance require the distribution of aerosol particle scattering cross- 
section n(o) [m _5 sr]. We assume a bi -modal, log-normal distribution of cross-section n(o) [David 
A. Bowdle, " n (o) Distributions," 1-31-90]: 


n(a) = 
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(27t) 1/2 /(X)oln[5 1 ] 
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where 



ln 2 8j 


ln 2 82 

/(X) = Oiexp 

2 

b, * 

+ xa 2 exp 

2 


is a normalizing parameter [m 2 sr 1 ], dj = 2x10 12 m 2 sr~ l , d 2 - 3x10 15 m 2 sr \ and 



(75) 


(76) 


where A x and A 2 are empirical coefficients [m -3 ] of the two aerosol modes, 8 t = 1.4 and 82 = 1.6 are 
the dimensionless geometric standard deviations of the two aerosol modes, and where n (o) - 
Bowdle’s n (a) [m -3 ] divided by a [m 2 sr~ l ] . For this example we chose x = 0.0003. 
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The SPM gives a biased estimate of P because the smallest aerosol particles (values of c ) pro- 
duce signal spikes below the threshold level, e<i„ 2 >. In this case ^ measured =0.185 p^ , i.e., there is 
a -81% bias in addition to the two error sources. The bias is independent of P for SPM if the shape 
of n(o) does not change, (i.e., a 1 ,02,8 1 ,5 2 , and x are constant). Note that there is no bias for VM 
measurements of p. The value of £ is a free parameter. Raising e will worsen the bias, cause less 
spikes to be seen in an observation time, but will reduce the number of noise spikes. Lowering £ 
will do the opposite. Moving to a tighter (closer) focus will lower (improve) the bias, but will 
increase alignment difficulty, and might cause the laser beam to change the properties of the aerosol 
particles. The effect this would have on <N spike > is not clear. 

Figure 10 shows the results for SPM. Both SPM noise error, and SPM random error assuming 
no bias due to the threshold level are plotted vs P . (The effects of bias due to threshold level on 
SPM random error could be calculated. A finite bias will increase SPM random error because the 
average number of spikes above the threshold decreases.) The average number of SPM signal events 
or spikes, > = N, that are greater than three times the average noise level (3 <t„ 2 >) in the 

measurement volume (V A T cos0 = 100m ) is shown in the figure. As long as the n (a) model is held 
constant, as in this case, N is proportional to p. Both sources of SPM error in p are proportional to 
P -1/2 . The curve N~ m , which is proportional to P -1/2 , is plotted for reference. In this example the 
random error is greater than the noise error and exceeds 10% when P < 3xl0“ 12 . As P approaches 
10 , the SPM number of spikes per observation time, N, approaches 1. It is very difficult to esti- 

mate P with a small number of signal spikes because of the large number of noise spikes. The noise 
spike correction of section 4.6 should be employed. Future analysis should include the effects of the 
noise spikes above the threshold. Increasing laser power will always improve SPM performance 
until the laser power begins to change the particles properties. The SPM analysis is only valid if 
indeed there are only single particle events. 

Comparing SPM to VM is a difficult undertaking fraught with many subtleties. Figure 1 1 is 
similar to Figure 10, but the N and N~ m curves are removed, and a VM noise error and two VM 
random error curves are added. The theory for VM error diverges whether one assumes operation in 
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a single-particle regime or in a many-particle regime (Gaussian statistics). The VM error curves in 
Figure 11 assume the many-particle VM regime except for the curve for VM random error in the 
SPM regime. Here, the many signal spikes below the SPM threshold level are mcluded in the VM 
estimate. This reduces the random fluctuations of the VM estimate compared to the SPM estimate 
which improves the VM accuracy. The random error for VM in the SPM regime may be better than 
the VM speckle regime. This can only be determined by calculating the boundary between the VM 
and SPM regimes. The VM random error in the speckle regime is independent of P and equals 
[T d !T] m . We assumed T d = signal decorrelation time = lps, and T - observation time = 1 s, 
yielding an error of 10 -3 (0.1%). The VM noise error is proportional to (T 1 and to B 1/2 [see Eq. 
(28)]. We assumed B VM = B SPM = 1.5 MHz for these calculations. In actual practice, one would try 
to reduce B VM as much as possible. The curve may be scaled for other values of B VM , e.g., if 
B vm = 150 kHz , the VM noise error is lower by 10 1/2 . 

For this particular example, SPM appears to be the best choice only for a narrow region near 
P= 10 r11 . For P > 3xKT u , the VM errors are lower, and for P < 3x10 l2 , the number of SPM 
spikes per observation time is less than 3. This eliminates SPM because the noise spikes would con- 
taminate the estimate. The SPM error is dominated by the random error component while the VM 
error is dominated by the additive noise component 

It must be remembered that the SPM measurement is biased low, while the VM measurement 
is unbiased; that the SPM and VM analyses are only valid in SPM and VM regimes, respectively; 
that the operational value of B VM may be much lower than B SPU ; and that systematic errors have 

not been addressed. 

The two methods of operation (VM and SPM) should be compared with separate optimization 
of parameters, and not necessarily equal focal settings or detection bandwidths, if the main goal is 
accurate measurement of P . However, if information about n(o) is also desired, then a tightly 
focused geometry is required, and the VM-SPM comparison should be done with the same system 
parameters. Without this theory it is not clear whether focused VM operation, which violates the 
Gaussian statistics assumption (many particles), is better or worse than collimated VM operation. 
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7. TRADE-OFFS 


7.1 VM Trade-offs 

The shot noise and random error (speckle statistics) contributions to total error for ideal VM 
operation are very small. The main sources of error are systematic. A focused beam has a factor of 2 
more system gain than a collimated beam. However, a focused beam senses a much smaller 
volume, may not have "speckle statistics," is more difficult to align, and may change or destroy the 
aerosol particles. Since systematic errors are likely dominant, a factor of 2 loss in gain with a col- 
limated beam could be worth the decrease in systematic errors. Improving beam quality at the 
expense of transmitted power may improve overall performance. 

12 SPM Trade-ofTs 

SPM operation is more complex and difficult to analyze. For useful operation, the system 
parameters must be chosen to provide single particle events over the specified range of {3. This 
requires a focused system where beam alignment may be difficult There are many trade-offs 
between beam geometry and system gain. The transmitted beam can be focused and the local oscil- 
lator beam collimated. This provides a larger range of SPM operation at the expense of system gain. 
A tightly focused system also requires higher bandwidth to observe the narrow brief pulses. This 

sacrifices SNR. Improving beam quality at the expense of transmitted power may improve overall 
performance. 
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8. FUTURE WORK 


8.1 VM 

1) Investigate temporal decorrelation time (we assumed 1 ). 

2) Investigate systematic errors: 

a) Beam geometry (collimated vs focused) 

b) Specifications on sensing volume and range 

c) Signal processing (rms detector, peak detector, power detector, linearity, background subtrac- 
tion, drift) 

3) ^random ) ^ single particlc behavior. Are there parameter regimes where VM estimates in 

the SPM regime have better accuracy compared to speckle statistics (many scattering particles) 
Le., collimated vs focused? 

4) Improved signal processing. What is the IF signal bandwidth (instantaneous and IF drift)? 

5) Is there an automatic gain control and what are the specifications? 

&2SPM 

1) Need reliable n (o) models. 

2) Calculate and A P( rflmto ^.2 for general parameters (beam size, focus, alignment, 

(J P 

beam offset). 

3) Calculate average number of spikes <N spike > for general system. 

4) QQrplatft distribution of spike amplitude and width. 

5) Determine SPM operating regions. 
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Figure 1. Geometry for a coherent detection laser radar system. 

An actual system would have overlap of the XMTR and 
BPLO beams at the target. 
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Figure 2. Geometry of single particle mode (SPM) operation 
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Simulation of IF power (the square of the IF voltage) 
with no noise for a moving point scatterer through a 
ussian transmitted beam. The average IF power (SPM 
lidar signal) is indicated by the dashed line. 
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9.0 


DEVELOPMENT OF A Ho : YAG SOLID-STATE PULSED LIDAR 

Spatially resolved, accurate measurements of atmospheric 
wind velocity are needed for many civilian, military, and 
commercial applications where the state-of-the atmosphere is 
important. Examples include weather forecasting [National 
Academy of Sciences (1980)], and the premonitory detection of low 
altitude wind shear by aircraft [R. Targ et al. (1990)]. This 
section documents the results of work performed during this 
effort toward the development of a pulsed, coherent, Tm,Ho.YAG, 

2 /xm laser radar system capable of making range-resolved 
measurement of atmospheric wind and aerosol backscatter profiles. 

Until recently, the mature CO 2 laser has been the almost 
exclusive type of transmitter employed in coherent laser radar 
(CLR) systems. The C0 2 laser and CLR technology has developed to 
the point where several CLR systems have been built for wind 
velocity and/or aerosol backscatter measurements. [see e.g., 
Hardesty, et al. (1988) and Menzies et al. (1984)]. The 
C0 2 -based LAWS instrument is currently being designed for 
space-based measurement of global wind fields [Beranek et al. 
(1989)]. Several issues relating to frequency stability, 
lifetime, pulse frequency chirp, and efficiency remain to be 
solved or demonstrated before C0 2 CLR transmitters can be space 


qualified [Hardesty et al. (1988)]. 

Solid-state laser technology is a rapidly advancing 
competing technology. With the advent of high efficiency, high 
power, semiconductor laser diodes as pump sources for solid-state 
lasers, all-solid-state CLR transmitters can be envisioned which 
have the advantages of low mass, small size, long shelf and 
operating lifetime, and the absence of consumables. In 
applications where atmospheric turbulence isn't too severe, the 
shorter wavelengths of solid-state lasers (0.7 - 3.0 fm) offer 
superior overall performance compared to the longer wavelengths 
of C0 2 lasers (9-11 /xm) . In addition, the continuous-tuning 
ranges of solid-state lasers are typically large enough to allow 
very low atmospheric extinction and/or multiwavelength 
differential absorption lidar (DIAL) measurements. 
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A 1.06-aot CLR system using Nd : YAG lasers was recently 
developed by Coherent Technologies, Inc. [Kavaya et al. (1989) ] . 
Velocity and range measurements from hard targets and atmospheric 
aerosol particles have been routinely demonstrated using this 
system at a remote site near Boulder, CO [Henderson et al. 
(1990A) ] . One of the chief deficiencies in the 1.06-/zm Nd : YAG 
technology is the lack of eye safety. 

The advances made in eyesafe 2 //m solid-state lasers during 
the mid-1980s led CTI to begin, under funding from the U.s. Air 
Force-Space Division, the development of the first eyesafe 
solid-state CLR system in 1987. Details of the 2-^m CLR 
development are included in the Air Force final report [Henderson 

et al. (1990B) ] . An overview of the 2-^m CLR is given in Section 
9.1. 

^ k^-"®^dboard master oscillator (MO) was developed under the 
Air Force funding and was utilized to demonstrate, for the first 
time, single-longitudinal-mode operation of a Tm,Ho: YAG laser 
[see Appendix E] . This breadboard was also used to "zero in" on 
a MO design which was functional. The NASA/MSFC funds provided 
during this effort were used to design, construct and test a 
"hardened" version of this MO which exhibits much better long and 
short-term frequency stability. This hardened MO is described in 
detail in Section 9.2. 

With the aid of funding from the Wright Research and 
Development Center (Wright-Patterson AFB) , the coherent 2-fim CLR 
system was recently relocated to the Table Mountain field test 
site located ~ 10 miles north of Boulder, CO. In Section 9.3, we 
describe preliminary field measurement results using the CLR. 

9.1 Overview of the 2-/zm CLR System 

This section contains an overview of the Tm , Ho : YAG coherent 
laser radar system. More detail about the master oscillator 
developed during this effort is provided in Section 9.2. The 
2-/m CLR system was primarily designed for measurements of 
atmospheric water vapor. 

A diagram illustrating the CTI Tm , Ho : YAG CLR system is shown 
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in Figure 12. The MO is a continuous-wave (CW) diode laser 
pumped single— longitudinal— mode (SLM) Tm,HotYAG laser. This 
oscillator, which we designed and built, utilizes a linear cavity 
with an intracavity Brewster plate for polarization control, and 
two solid fused silica etalons for frequency control. With the 
MO laser crystal thermo-electrically cooled to ~ 230 K, the SLM 
output power is 50 mW near 2090 and 2097 nm. The wavelength 
can be tuned continuously between 2086 and 2100 nm. By changing 
the reflectivity of the output coupler, tuning near 2121 nm has 
also been achieved. The linewidth and short-term frequency 
jitter (seconds) of the MO is < 1 MHz, and the long-term 
frequency drift (~ 1 day) is ** 300 MHz peak— to— peak . 

The output of the MO is first passed through a Faraday 
isolator which isolates the MO from the rest of the system. 
Additional isolation from the slave oscillator is provided by a 
98% reflector. A portion of the reflected power serves as the 
local oscillator (LO) beam. The MO power that is transmitted by 
the 98% reflector is frequency shifted by 27 MHz using an 
acousto-optic frequency shifter (AOM in Figure 12). The 
frequency— shifted beam is then mode— matched and injected into the 
Q-switched slave oscillator. The f lashlamp-pumped Cr , Tm , Ho : YAG 
laser crystal in the slave oscillator is water cooled and 
operates at approximately room temperature. The injection-seeded 
slave oscillator is capable of producing 50 mJ, ~ 150 ns, SLM, 
TEMo o pulses at 3 Hz . The required frequency match between the 
slave oscillator and the MO is maintained by an automatic servo 
system which controls the slave cavity length using an end mirror 
mounted on a piezo-electric translator (PZT) . More detail about 
the performance of the injection-seeded transmitter is provided 
in Appendix F. 

The SLM output of the slave oscillator is incident upon an 
eccentric-pupil Dall-Kirkham telescope which expands the beam to 
~ 20 cm and transmits the pulse into the atmosphere. 
Backscattered radiation from the atmosphere is collected by the 
telescope and mixed with the LO beam using a beam splitter whose 
reflectivity is ~ 15%. The coherently-mixed signal and LO beams 
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are then focused onto a 100 /im-diameter room- temperature 
extended-wavelength InGaAs detector. The quantum efficiency of 
this detector at 2.1 jm was measured to be 0.69 electrons/photon. 
The heterodyne signal from the detector, which is in the 
bandwidth between * 5 and 50 MHz (target velocities from ~ 23 m/s 
to ^ —24 m/s), is amplified, filtered, then digitized every 10 ns 
using one channel of a LeCroy 9400 dual-channel digital storage 
oscilloscope (DSO) . The digital data is then transferred to a 
DEC MicroVax II digital computer using a general purpose 
interface bus (GPIB) . The raw data is then processed using very 
flexible processing software, to extract Doppler velocity (first 
moment) , velocity width (second moment) , and backscatter (zero 
moment) information. 

The CLR is constructed on a 4' x 8 7 Newport optical table. 
The laser transmitter consisting of the MO, Faraday isolator, 
AOM, and slave oscillator, is mounted on a separate 1 ' x 5' 
breadboard. Figure 13 consists of two photographs showing 
different views of the complete laser transmitter. The scale in 
the photos is given by the 1” hole pattern of the 1' x 5' 
breadboard . 

9.2 Frequency Stabilized Master Oscillator 

As described in Section 9.0, a breadboard oscillator was 
used to demonstrate that SLM oscillation could be achieved. This 
hardware allowed us to try many different resonator 
configurations to arrive at an operating condition which 
performed best. The performance we were able to achieve with 
this breadboard hardware is described in Appendix E and in the 
final report to the Air Force [Henderson et al. (1990B)). A 
photograph of the breadboard hardware is shown on the top half of 
Figure 14. 

The performance data collected with the breadboard-level 
oscillator was used to design and construct a highly stabilized 
SLM laser (Figure 14 bottom) . In addition to mechanical and 
thermal improvements, a significant change in the basic cavity 
design (as compared to the design presented in Appendix E) was 
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incorporated as a result of earlier tests. The resonator 
configuration chosen for the stabilized laser is shown in Figure 
15. The intracavity lens Li allows for a "collimated" leg inside 
the cavity between the lens and the output coupler M 2 , where the 
Brewster plate polarizer P and tilted etalons Ei and E 2 could be 
installed. This effectively reduces the walkoff loss experienced 
in those elements, which reduces overall circulating intracavity 
losses and permits higher tilt angles of the etalons before 
excessive loss is incurred. This in turn assures a very clean 
TEMoo spatial profile in the output beam, as all intracavity 
optics (and their inevitable surface reflections) are capable of 
being tilted well off of the beam axis. This technique also 
reduces output beam astigmatism due to the tilted Brewster 
polarizer to a negligible level. Very high spatial beam quality 
in the LO source is essential for maximum heterodyne efficiency, 
and could not be compromised in the laser design. Details of the 
resonator optics follow. 

The flat resonator mirror Mi is a high reflector at 
wavelengths near 2100 nm and transmits ~ 88% of the diode laser 
pump light near 785 nm. Resonator mirror M 2 is also flat and 
serves as the output coupler, transmitting ** 2.5% of the laser 
light near 2100 nm. The cavity mirrors Mi and M 2 are separated 
by " 117 mm. The intracavity lens Li has a focal length of ~ 75 
mm, is located ~ 75 mm from Mi , and is AR-coated at 2100 nm. The 
uncoated solid-fused silica etalon E 2 is 0.09 mm thick. Ei is a 
0.5 mm thick solid-fused silica etalon with the surfaces coated 
for 18% power reflectivity near 2100 nm. The Tm:Ho:YAG laser rod 
is 4 mm in diameter and 4 mm long. The end faces of the crystal 
are polished flat and AR-coated for wavelength near 2100 nm. The 
AR-coating reflects < 0.5% of the 785 nm pump light. The laser 
rod is cooled to ~ -40° C using a 3-stage thermo-electric cooler. 

A great deal of effort was expended on the opto-mechanical, 
thermal, and acoustical engineering issues associated with 
producing a very frequency-stable oscillator. To assure absolute 
frequency jitter of less than 1 MHz out of 1.4 x 10 1 4 Hz (i.e., 
2.1 ^m wavelength), requires dimensional stability of the 
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JV 


resonator mirror spacing on the order of 
length of ~ i 17 ^ . In addition. 


w 8 i (for the cavity 
it is desirable to limit 


long-term frequency drift due to thermal and/or pressure changes 
in the environment to a very low level. Negligible frequency 
drift must be achieved in the MO/LO both to assure that the 


effects of atmospheric absorption remain constant during laser 


radar operation, and to simplify the frequency-locking 
requirements in the injection— seeding subsystem. 


For these reasons, we designed and built a box-like 
resonator structure made from super-invar, an ultra-low thermal 
expansion steel alloy, shown in Figure 14 (bottom) . Both 
resonator mirrors are located against solid shoulders milled into 
either end of the common structure, assuring immunity against 
mechanical vibration and, to first order, thermal changes in the 
environment. The entire super-invar body is temperature 
stabilized using a small closed-loop constant-temperature water 
circulator. This arrangement also permits cooling and thermal 
stabilization of the thermoelectric (TE) heat pump module used to 
cool the laser rod. This TE cooler module is powered by a very 
low noise (a few mV p-p) DC power supply, further assuring that 
thermally-induced refractive index fluctuations in the laser rod 
are held to a minimum. 


To provide isolation from acoustic vibrations and pressure 
changes in the environment, the entire resonator body structure 
is hermetically— sealed and filled with dry nitrogen. The 
n ^ ro 9 en backfill also prevents frost buildup on the YAG rod, 
which is typically operated at -40° c. The stabilized laser's 
output power and frequency tuning, as observed with the 
quarter-meter monochromator and complementary spectrum analyzers, 
produced data very similar to that described in Appendix E for 
the breadboard- level oscillator. 


As expected, the engineering refinements described above 
have resulted in very low short-term jitter and drift in the 
laser output frequency. By optimizing the finesse of a scanning 
confocal spectrum analyzer (free spectral range of 300 MHz) , it 
was possible to resolve and characterize sub-megahertz frequency 
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fluctuations in the stabilized laser output. Figure 16 shows 
typical frequency jitter and drift magnitudes for one-second and 
five-second operating periods. The horizontal dimension of the 
oscilloscope traces represents 1.13 MHz per major division (200 
^s) ; frequency jitter is observed as a horizontal movement of the 
interferometer transmission profile over the scope camera 
exposure time. The width of the transmission profile is due to 
the resolution of the 300 MHz interferometer. Typical jitter and 
drift values observed were less than 0 . 5 MHz peak-to-peak for up 
to 5 second periods, and correspondingly much less for the 
short-term (10's of /zsec) round trip times of interest for 
atmospheric coherent laser radar. Long-term frequency drift has 
been found to be * 300 MHz /day peak-to-peak in this laser, 
resulting in extremely stable, mode-hop-free SLM operation. It 
should also be noted that jitter and drift recorded in these 
photographs are contributed to further by the scanning 
interferometer diagnostic itself. Finally, the photographs were 
taken with the laser and interferometer mounted on a standard 
optical table which was not vibration isolated. The optical 
table was located inside a relatively noisy laboratory on the 
third floor of an office building. 

9.3 Preliminary Field Measurement Results 

The system was very recently moved and integrated at our 
remote field test site ~ 10 miles north of Boulder, CO (Building 
T-l at Table Mountain - a DOC facility) . Only a small amount of 
CLR data has been taken at this time. We include here examples 
showing some of the preliminary results which exhibit the 
fundamental operation of the CLR system. 

For all of the examples shown below, the transmitted pulse 
energy of the CLR was ~ 20 mJ. The CLR was running at a PRF of 
3.2 Hz. A LeCroy 9400 DSO was used to capture the temporal data 
shown in the following figures. Figure 17 shows the temporal 
profile of the 20 mJ pulses. Note that the full-width at 
half-maximum (FWHM) duration of the SLM pulses is ~ 220 ns. 
Figure 18 shows the heterodyne signal due to the single-shot CLR 
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return from a mountainside located * 14.9 km down range. The 
telescope was focused at infinity for this data taken on October 
4, 1990. Figure 19 shows an example of a return from clouds 
located ~ 13.8 km from the CLR taken on October 4, 1990. When 
this single-shot data was taken, the telescope was focused at 
2 km. 

Figure 20 plots the power spectrum of the return from a 
cloud located approximately 10 km from the CLR. A 256-pt. 
fast-Fourier-transform (FFT) centered at 9.89 km was used to 
compute the power spectrum. Note that the radial velocity of the 
cloud is ~ 7 m/s toward the laser radar. Figure 21 gives an 
example of a wind measurement made on October 7, 1990. The 
telescope was focused at 2 km and the beam was aimed to the north 
at ~ +8° from horizontal. Figure 21 is the plot of the power 
spectrum of aerosol return in a thin slab of atmosphere centered 
at 3.75 km. A 256— pt. FFT was used to generate the spectrum. 
Note that the radial wind velocity is ~ 8 m/s toward the CLR. 

To date, atmospheric wind velocity has successfully been 
measured to horizontal ranges exceeding 20 km on hazy days and 5 
km on very clear days. When the lidar beam is pointed 
vertically, the highest altitudes above ground level for which 
aerosol return has been observed to date is ~ 4 km. Figure 22 
shows the A-scope display of aerosol return signal for 111 shots 
averaged. For this data taken on 24 Oct. 1990, the lidar was 
aimed approximately horizontal and focused at m. The receiver 
bandwidth was ~ 41 MHz. The A-scope display is the result of 
squaring the IF aerosol signal from each shot and averaging 111 
consecutive shots. The lower trace of Figure 22 defines the 
background noise and is the result of averaging 111 shots with 
the telescope output blocked. When narrowband FFT processing is 
utilized, the signal-to-noise levels shown in Figure 22 are 
sufficient for wind measurements beyond 20 km (> 13 0 /zs) . 
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Figure 13. Injection-seeded Cr , Tm , Ho : YAG transmitter, comprised 
of master oscillator, slave oscillator, and 27 MHz 
AO modulator. 
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Figure 14. CW Master Oscillator. Top: Breadboard-level 

hardware configuration. Bottom: Final stabilized 
oscillator with top removed from resonator structure. 
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5 seconds 



Figure 16. Typical frequency jitter and drift measured during 
1 second and 5 second observation times using a 
300 MHz free spectral range interferometer. The 
combined jitter and drift is <0.45 MHz peak to peak. 
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Figure 17. Intensity vs time of the transmitted 2.1 um 
CLR pulses with the pulse energy set to 
~ 20 mJ. 
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Figure 18. Single-shot CLR return on 4 Oct. 1990 from 
a mountainside located 14.9 km down range. 
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Figure 19. Single-shot CLR return on 4 Oct. 1990 from 
clouds located at ~13.8 km range. 
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Figure 20 . Power spectrum of a single-shot cloud return 
located ~9.9 km from the CLR. 


57 



s 


m 

5 

111 


i 


12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 


POWER SPECTRUM AT RANGE 3.750 KM 

(AWAY FROM UDAR) VELOCITY (M/S) (TOWARD UDAR) 

25 20 15 10 5 0 -5 -10 -15 -20 



■ 


• 

- 

• 

— — — — — - 

- 

• 

- 

• 

IBBBH 

- 

BflBBE 

• 

IBBB^fll 


BIBBB 

■BEEBES 

BEEBBI 


BBBB 




FREQUENCY (MHZ) 


Figure 21. Power spectrum of aerosol return from 3.75 km 
range taken on 7 Oct. 1990. 
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Figure 22. Upper Trace: A-scope display of aerosol return 

from 0-30 km for 111 shots averaged. Lower 
Trace: Background noise level. 
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10. SUMMARY AND CONCLUSIONS 

VM is much simpler and more promising than SPM. The SPM 
analysis is much harder than the VM analysis. The non-systematic 
errors for VM (shot noise and speckle contributions) are less 
than 1% for 0 > lO'^m^sr' 1 with a 10 s measurement time. 
Therefore, VM systematic errors will dominate. Priority should 
be given to perfecting VM operation, e.g., ensuring small 
systematic errors by optimal choice of system parameters 
(sensitivity analysis), improving alignment methods, improving 
signal processing methods, improving field calibration 
techniques, and improving beam quality. A crucial element of VM 
and SPM operation is accurate conversion of IF signal to IF 
power. Alternativces to a SAW processor should be investigated 
and may prove to be cheaper, more accurate (linear) , and more 
reliable (less drift, etc.). 

A Tm, Ho : YAG 2-/*m laser radar system has been constructed and 
preliminary measurements demonstrating the first laser radar 
returns using an eyesafe solid-state CLR have been made. 
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APPENDIX A. Volume Mode Estimation of f) 


The average coherent lidar power <i/>[A ] is given by 

1T1 ,^ 2 - - 

<//> = 2 


T\G D e 


hv 


j j M s <? U %0)W( *i)W' ^2 )Elo * C^i.O )E lo (V 2 ,0 )d^ 2 (Al) 


t] = detector AC quantum efficiency [electrons/photon] 

G d = amplifier gain [dimensionless] 
h = 6.626xl( T 34 [Js ] = Plancks constant 

v = = frequency of the optical field [s _1 ] 

Euo (V,0) = LO field back-propagated to receiver plane [Wm" 2 ] 1/2 
W = transfer function of receiver lens [dimensionless] 

(V^O) = <Es (?uO)E s * (V 2 ,0)> = mutual coherence function of received field at the receiver 
[Wm~ 2 ] 


For randomly distributed aerosol particles, the mutual coherence function of the total backscattered 
field is the sum of the mutual coherence functions from each aerosol particle. The mutual coherence 
of the backscattered field from an aerosol at (ftjt ) with cross section a [m 2 sr~ x ] is 

Ms (ft* )G 0W )G * CW )> (A2) 

where G (V#,/? ) is the Green’s function [m -4 ]; i.e., the field at (V,0) due to a unit point source at 
(pji ). The mutual coherence function of the total field is the sum of the mutual coherence functions 
of all the scattering aerosol particles, i.e., 

M s (Vi,V 2 ,0) = J J f\ 2 Gn (Gjt )<Ij (pjt )G )G * (V 2 ,jt/? )>d GdftdR (A3) 

0--0 

where /t(o^?)[m -5 sr] is the number of aerosol particles per unit volume per unit o at range R. If 
n (Gjt ) is independent of range R, then 
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where 


(3 = jan(a)da 
o 


is the aerosol backscatter coefficient [m~ l sr 1 ]. The average coherent lidar power is then 

"tiG € ' ** ** ** 

<i/> ■ 2p ~j~~ f J J IfchtfJnGCtajClG'OfrfJtfr 

k J v - - - 


C?2 )Elo * <yi,0)£ LO 


It may be written 


<t’s 2 > = fiGyP t 


where Gy is the system gain [A 2 mW l sr], i.e.. 


G v = jG v (R)dR 
o 


-r-j j J <I T (tJt)0(^„TtJ»0’ 

n V ij ; w 

k J 

VWi)W* * (^i.0)£lo 

is the coherent lidar system gain density [A 2 ^rW~ J ]. 
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APPENDIX B. Single Particle Mode Estimation of p 


The mutual coherence function of the backscattered field from a single particle of cross section 
a [m 2 sr _1 ], which is far enough away that the received field approximates the field from a point 
source under the parabolic approximation, is given by 


m s ew»- 


o/Cftfi) 


exp 




(P-W 


(Bl) 


where 

j}=[p x j)y] = transverse trajectory of the particle [m] at range R, 


/ r aW = 


(2xRr — — 


* 

J J E T (p h 0)E T * 


dtt x d^ 2 > (B2) 


If = [«x ] is the transverse coordinate [m] at the transmitter plane perpendicular to the R axis, 

E r (if,0) = transmitted field \Wm ~ 2 ] 1/2 at R = 0, ~~ 

/ r (tf,0) = E T (tt,0)E T * (Vt,0) = transmitted irradiance [Wm~ 2 ]. 


The maximum lidar signal S, [A 2 ] for the jth particle event will occur for some coordinate 
(Px Pym J * ) °f the trajectory, where p x is the displacement from the laser beam axis. This maximum 

signal is given by 


Sj(PxfymJt) = 


2P t g 

T[Gpe 

R 2 

hv 

* J 


2 M - 




WftW* (?2 )Elo * (tiJMiO 


(B3) 


where Jf m = (p x p ym ) is the transverse coordinate of the m a xim u m signal for the jth particle (see Fig 
ure 2) and J T (pJl) = h (ft* ) /p T fa -2 ] is * c normalized irradiance. 

The proper statistic for estimating p is the sum S [A 2 ] of the maximum (or peak) signals: 


N 


j = 1 


(B4) 
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For large N (many particle events), this statistic will be very close to its ensemble average <S>. 
(Since is assumed constant over the observation volume, the statistical process is ergodic and sta- 
tionary.) The observation volume is determined by the observation time T[s], During this time the 
lidar system will have traversed a distance V A T through the atmosphere. We obtain <S> by per- 
forming an ensemble average over all the aerosol particles contained in the observation volume: 


QO mm V^Too m oo V A T<m om mm 

< S > = J J j jSj(p x j) Ji)n(o)dadydp x dz =2 P T j f j f f f 

0-00 o - 0 0-- 


T]G D e 


^2 


hv 


an (a) 


{t ®rn£) 

R 2 


exp 


-Vi> 2 -(ft. -V] WpJW'ctJMu, dljdnfydp,* (B5) 


The integrand in Eq. (B5) is independent of y so the y integral may be evaluated, producing V A T. If 
the laser’s optic axis is at angle 0 to the z axis, then z = R cos9 . Thejingle hg avi *; nnrT 
the particles trajectories will be 90® - 0 ._Ths integration OV€rz~*is more conveniently performed by 
chan^ngj^i^le^tismf^r= dR cos0 and integrating over R . The a integration may be performed 
if we assume n (a) is constant over the dominant region of integration. Then 


<5 > = p P T G s V A T cos0, (B6) 

where G s [A 2 W~ l sr ] is the total system gain given by 

« M 

Gs = J J G s (p x Jl)dp x dR. (B7) 

o -- 


(The definition of p in Eq. (B6) is the same for direct detection and heterodyne detection volume 
mode measurements.) G s (p z Ji) is the system’s transverse area gain density [A 2 W~ l m~ 2 sr ] for a 
scattering particle at (pjt ) given by 


G s (p xt R) = 2 


r\G D e 


hv 
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J J 




R‘ 


-exp 


^-[(Hn-^l) 2 -0^,-^2) 2 ] 


W&0W' * (1^,0 )E W (^,0 


(B8) 
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Section 4 assumed that the calibration of the lidar system utilized apparatus that shot individual 
tesT particles^hrough the beam with known test parameters a c J*c fxc .and Re- Both p x and R 
were assumed adjustable over a sufficient volume. This appendix treats the case when the particle’s 
value of p K is not known, but the probability density function (PDF) f (p^) of the value p x is uni- 
form. Specifically, 


1 

if -L/2 < p x < LI2 


f(Pxc) 


=< 


0 


otherwise 


(Cl) 


We define the Monte Carlo statistic Z to be 


Z(Pxc) = 


L 

P C <*C N C 


N c 

Z^c(Pxc»^c) 


7=1 


(C2) 


where S c (p x &c ) is * c maximum lidar signal power for the jth calibration particle event, p^ is the 
random displacement from the laser beam axis, and N c is the (large) number of observed events. 


The average of the statistic Z is 


L/2 N c 


<Z(p K )> = \ f(Pxc)Z(Pxc )dPxc = p c a c N c J rl 


(C3) 


If we assume L is much larger than the beam width at Rq , the integration may be considered from 
. oo to <x>, and, referring to Eqs. (64) and (65), 

<Z> = G S (R c )■ <C4) 


Equation (C4) is used to obtain the system gain G s in Eq. (36) to calibrate (3. Any deviation from the 
uniform PDF behavior of p K will produce an error in (3. 
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ABSTRACT 

The SNR equation of a CW, monostatic, coherent laser radar 
system is examined for the case of a distributed aerosol target. 
Calculations and plots are presented showing the location in 
range and the volume of the laser radar measurement. 

KEY WORDS 
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1. Introduction 

Monostatic, coherent detection laser radar (LADAR or LIDAR, 
light detection and ranging) systems with continuous-wave (CW) 
laser sources, with the atmosphere serving as a distributed 
target, have proven to be important remote sensing tools for many 
applications such as measurements of radial or line-of-sight 
(LOS) wind velocity, 1-6 true airspeed of aircraft, 7-8 transverse 
wind velocity, 9 aircraft trailing vortices, 10 and atmospheric 
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aerosol backscatter coefficients . 11-14 Comparisons of CW laser 
radar wind measurements with other wind velocity monitors have 
been performed, 3,6,7,9,15,16 as have investigations of intensity 
calibration issues .11*12,17,18 

A separate but very important consideration when discussing 
CW coherent laser radar (CLR) measurements based on backscatter 
from atmospheric aerosol particles is the question: "Exactly 
where along the optic (range) axis is the measurement being 
made?" This question comprises the issues of 1 ) whether the 
range resolution of the measurement is too small or too big, 2 ) 
whether a high backscatter region far from the laser radar (e.g., 
cloud or ground) corrupts the desired backscatter or wind 
velocity measurement , 18 3 ) whether a region of large wind shear 
far from the laser radar corrupts a wind velocity measurement by 
broadening the signal spectrum, and 4 ) over what ranges must a 
hard calibration target be used to achieve a desired 
characterization of the laser radar. The laser beam diameter vs 
range is important to ascertain eye safety and whether the beam 
irradiance may alter the properties of the aerosol particles. 
Also, the laser radar measurement volume is important since too 
small a volume may violate the statistical assumptions (many 
particles) of the laser radar theory, or may cause a measurement 
bias by missing scarce yet important large aerosol particles that 
contribute to the backscatter coefficient. The measurement 
volume is also important in explaining differences between CW 
laser radar measurements, and measurements made with other 
instruments. To our knowledge, published discussions of these 
effects have been limited to the range interval yielding, e.g., 
50 % of the signal under tight focusing conditions . 1 ' 3 ' 7 ' 10-12 ' 16 

In this paper we solve the laser radar SNR equation for any 
percentile range of the received signal, and therefore for any 
definition of range resolution; and we also solve for the 
measurement volume of a CW coherent laser radar system under any 
focusing conditions. 
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2 . 


SNR Equation 

We start with the general untruncated Gaussian beam 
monostatic, pulsed, shot-noise-limited CLR signal-to-noise ratio 
(SNR) equation and assume: 1) CW transmit laser power P [W] , 
2) negligible atmospheric extinction, 3) negligible atmospheric 
refractive turbulence effects, 4) constant aerosol backscatter 
coefficient 0 [m _1 sr -1 ] vs range R [m] , 5) matched diameter and 
phase curvature F [m] of the untruncated Gaussian transmitted and 
back propagated local oscillator (BPLO) optical fields, 6) 
negligible beam truncation by the transmit/receive lens or 
telescope, and 7) a photovoltaic detector. Under these 
conditions the dimensionless SNR equation becomes 1 9 


SNR 


^cw£ 

hz/B 


r <D 
- 0 



dR 

a -S ) 2 * 2d c 


4A‘ 


” P cw^ 

hz/B 



¥ R 

R 2 + (1 - f) 2 R^ 


( 1 ) 


where r) [electrons/photon] is the detector quantum efficiency at 
the signal frequency, h = 6.626 10‘ 34 [Js] is Planck's constant, 
v [Hz] is the optical frequency, B [Hz] is the receiver 
bandwidth, A = cfv [m] is the wavelength of the optical field, c 
[m/s] is the speed of light, and D Q [m] is the e" 2 (14%) 
intensity diameter of the transmitted and BPLO Gaussian fields. 
The Rayleigh range [m] , 



(2) 


is defined for Gaussian beams as the distance from the beam waist 
(minimum diameter ■ D q ) to the point where the beam area has 
doubled. 20 Alternatively, for light diffracted by an aperture of 
diameter D q , ranges much larger than Rj^ will satisfy the 
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Fraunhofer or far-field assumption, and the field at range R will 
be the Fourier transform of the field in the aperture, 21 changing 
only in size but not in shape with increasing R. 

The two terms in the denominator of Eq. (1) reduce the 
coherent detection SNR. The first term represents the effect of 
the finite receiver solid angle, and the second term represents 
both the speckle effects from the size of the illuminated diffuse 
aerosol target, and the phase-front mismatch (PFM) at the 
receiver when R ^ F. The integrand I in Eq. (1) is the 
"weighting function" of the received signal along the range axis 
under the assumed conditions. For small R (R « F and R « R^) , 
the integrand becomes independent of R and equal to R^ -1 . As R 
decreases, the gain in SNR, due to increasing receiver solid 
angle, exactly cancels the SNR loss due to speckle/PFM. For 
large R (R » F cje R » R^) the integrand decreases as R -2 due to 
the decreasing receiver solid angle. 

For R - F, the integrand equals 


I(R = F) 




( 3 ) 


which is the proportional to the effective solid angle of the 
receiver. 

The range to the peak Rp [m] of the SNR weighting function 
is found by setting the derivative of the integrand in Eq. (1) 
with respect to range R to zero, yielding: 



F 

1 + F 


2 

N 


( 4 ) 


where 


F 


N 


A 


FVRr 


( 5 ) 
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is the focal range normalized by the Rayleigh range 

(dimensionless) . For small F N (F << *R>' 

Rp « F, (6) 

and the peak SNR occurs at the focal range (i.e., tight 
focusing) . For large F N (F » *R>' 



(7) 


and the peak SNR range is proportional to F" 1 , i.e., the peak 
originates close to the laser radar. One result of this is that 
for laser radars focused far . away (F » Rp) , the very term 
"focus" loses its meaning. The range to the peak SNR goes to 
zero for both very small and very large F. Since Eq. (4) is also 
the range to the Gaussian beam waist, 20 we conclude that the peak 
SNR always comes from the beam waist position. 

For R * Rp, the peak SNR will be proportional to the 
integrand in Eq. (1) which becomes 


I(R - Rp) 



( 8 ) 


For small F N (F « Rp) , I(R - Rp) « RpF -2 and the P eak SNR 9 ets 
very large since receiver solid angle SNR losses and speckle/PFM 
losses are simultaneously low. For large (F » Rp) , I(R = Rp) 
« r^ 1 , a constant value identical to the small R value of I, 
since R p is small and we assumed R = Rp. The ratio of the SNR at 
the focal range to the peak SNR is found by combining Eqs. (3) 
and ( 8 ) : 


SNR(R=F) = **R = 1 

SNR ( R=R p ) r 2 + F 2 x + 
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For small F N (F « Rp) , the ratio approaches unity, indicating 
that the peak SNR occurs at the focal range. For large F N (F » 
r^) , the ratio approaches zero as F~ 2 , indicating that a 
negligible contribution to SNR comes from the focal range. 

The value of F which maximizes Rp is found by setting the 
derivative of Eg. (4) with respect to F to zero, yielding; 

Fp - Rp (10) 

The value of R when F - F is found by substituting Eq. (10) 

p P 

into Eq. (4) : 


vv 


*r 


(ii) 


When F = Rp, the peak SNR (and the beam waist) is at its greatest 
range and this range is Rp/2. Focus settings smaller or larger 
will cause the peak SNR to move closer to the laser radar. As F 
is increased from 0 towards ®, the peak SNR decreases and moves 
to larger ranges until F = Rp. Then the peak SNR continues to 
decrease and moves back to smaller ranges. For small F, the 
illuminated aerosol spot at R * F is very small, and the losses 
due to finite receiver solid angle and speckle/ PFM are both 
small. On the other hand, the losses due to speckle/PFM increase 
quickly away from the focus . Therefore the SNR is high and peaks 
sharply at Rp ~ F. For large F, the effects of speckle/PFM loss 
compete with the effect of R 2 , producing a broad SNR curve with a 
peak collocated with the beam waist at Rp < F. When the focal 
range is set to the Rayleigh range, the integrand value at the 

Rayleigh range is I (R - F - Rp) = V * The integrand takes on 
this same value for very small R, from which we draw the 
remarkable conclusion that, e.g., the first cm of range 
contributes the same SNR as a 1 cm slice of range at R “ 
which may be many 1cm away. Put another way, the solid angle SNR 
loss at R = Rp exactly equals the speckle/PFM SNR loss at small 

R. 
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As shown by Sonnenschein and Horrigan, 22 the integral in Eq. 
(1) may be solved analytically. Eq. (1) may be written as 


SNR 


l P cw^ 



1 







( 12 ) 


where 


Rn * R/Rr ( 13 ) 

is the range normalized by the Rayleigh range (dimensionless) . 

Letting a = 1, b = -2/F N , c * 1 + 1/F 2 , and A = 4ac-b 2 = 4, we 

may look up the integral 2 3 

f — 5— - 2A’ 1/2 tan” 1 !" (b + 2cx)A” 1/2 "| ( 14 ) 

J a + bx + cx* L J 

for A > 0. Using Eq. ( 14 ), Eq. (12) becomes 

SN * (•» - IHs** i 1 + 4 - 1”-' 1 (47]] < 15 > 

where the argument oo for SNR indicates that all ranges from 0 to 
gd are included. 

Note that only the second term in the brackets contains any 
dependence on D Q or F, the beam diameter and focal range. For a 
sufficiently tight focus, F N « 1, the expression in the brackets 
equals 2. For Fjj >> l, i.e., a collimated beam, the expression 
in the brackets equals 1 . Therefore, the beam diameter and focal 
range may affect the SNR only within a factor of 2 (3 dB) . A 

tightly focused CW laser radar will have twice the SNR of a 
collimated CW laser radar in a uniform atmosphere. (Recall that 
we are neglecting atmospheric extinction and refractive 
turbulence effects.) 
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3 . 


Comparison to Past Work 

Eq. (15) may be directly compared to Eq. (24) of 

Sonnenschein and Horrigan 22 . The RHS of their Eq. (24) must be 
multiplied by 2 due to an unreported error. Also their variable 
"R” is the e' 2 intensity radius and corresponds to our D/2. 
Finally, our Eq. (1) is for a photovoltaic detector as is their 
Eq. (24) when multiplied by 2. With these caveats, our Eq. (15) 
equals their Eq. (24) . 

4. Wavelength Dependence of SNR 

Since v = c/A, Eq. (15) contains an explicit A 2 dependence. 

The detection bandwidth B will also be proportional to A -1 if 

allowance is made for a fixed atmospheric velocity width Av 

[m/s], or for a maximum unknown radial velocity of the aerosol 

particles, V [m/s], either toward or away from the laser 
max 

radar. Ignoring any A dependence of r} or P cw # and neglecting the 
terms in the brackets, we see that SNR(a>) a /3X 3 . The A 
dependence of varies greatly in nature but lies between A 0 (Mie 

limit) and A -4 (Rayleigh limit). A commonly used approximation 
that 0 a A -2 would result in SNR(®) a A, favoring larger 
wavelengths. Of course, larger values of A may lower SNR(m) by 
as much as 3 dB through the bracketed terms in Eq. (15) . 

(Strictly, the above reasoning applies only to aerosol 

backscatter measurements. Wind velocity error expressions 
include both SNR and other terms that also depend on A . ) 


5 . Cumulative SNR 

The dimensionless cumulative SNR is defined as 


SNR 


CUMSNR 




( R N ) 


SNR(o) 


(16) 


where SNR(a>) is given in Eq. (15). Using Eq. (12) 
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SNR 
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* p cw^ A ^ 


hi'B 
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dz 
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1 + 


N 


(17) 


Employing Eq. (14) to solve Eq. (17) produces 

SNR <V - %^{ tan ' 1 [- -j; + [ l + pH] + tan '%]} 

" N 

We now use the trigonometric identity 23 


(18) 


tan” 1 x + tan 1 y = a + tan 1 i - xy ~] 


(19) 


where o * 0[rad] if xy < 1, ct = x if x > 0 and xy > 1, and a = -x 
if x < 0 and xy > 1. We set 


and 


x = 


N 


y(V 


i + 


N 


*N - 


N 


( 20 ) 


( 21 ) 


yielding 


xyfRjj) 


1 



m ^ 

i + 1 

F 

N 

V. J 



1 + 2 
F Z 
N J 


*N - 


N 


( 22 ) 


In the limit of a tight focus, F N << 1, xy(Rjj) ~ RjjF^ 3 which will 
be greater than 1 except for very small values of Rjj. In the 
collimated limit, F N » 1, xyfR^ ~ R^ 1 which will be less than 
1 except when R > F. Using Eqs. (18) and (19) we find 

SNR <V " hX 
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( 23 ) 


where a is shown as a function of Rjj since a depends on the value 
of xy^). Combining Eqs. (15), (16), and (23), we obtain 


CUMSNRfRjj) 




+ tan 


-1 



*N 


1 - 

‘V 



Lf n J 

- 


i r N j 


f + tan 


(24) 


CUMSNR(Rjj) should range from 0 to 1 as Rjj goes from 0 to ®. For 

very small xy(R N ) ~ , a - 0, and CUHSNR(O) - 0. For very 

large R^j, a ■ t, and using 

tan -1 (-z) = - f + tan” 1 ^) (25) 


for z > 0, we find CUMSNR (®) * 1, as expected. 

Figure 1 shows CUMSNR vs normalized range R^ for normalized 
focal ranges F N of 0.1, 0.5, 1, 5, and 10. The range at which 
CUMSNR * 0.5, the 50th percentile range Rjj( 50%) , is always 
smaller than F N , with the difference increasing with increasing 
F N . (Inclusion of atmospheric extinction and refractive 
turbulence in the theory would increase the difference between 
R n ( 50%) and F^. ) The range resolution, e.g., the range interval 
over which cumulative SNR increases from 5% to 95%, increases 
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with increasing values of F N> Note that for some values of Rjj, 
the cumulative SNR does not monotonically decrease with 
increasing values of F^. The curves may be translated into 
actual range R and focal range F values through knowledge of the 
Rayleigh range of a particular laser radar system. 

6. Percentile Ranges and Range Resolution 

As stated in the introduction, we would like to know the 
range to any percentile of the SNR. For a desired value of 
CUMSNRfRjj) , e.g., 5%, 25%, 50%, 75%, 95%, etc., Eq. (24) may be 
inverted to solve for R^: 

F tan? 

R N (CUMSNR) = F N + tan? (26) 

where the angle y [rad] is given by 

7 « (CUMSNR) j^T - tan _1 F N j (27) 

Eqs. (26) and (27) are an analytic solution for RjjfCUMSNR) . 

Figure 2 shows five normalized percentile ranges R^j plotted 
vs normalized focal range F^. The R JJ (95%) curve (dashed line) is 

shown divided by 10 to fit on the plot. Note the spread of 

curves once the focal range F exceeds the Rayleigh range Rj^ 
(i.e., logio F N > 0) . Rj^ becomes equal to F N for small F^. For 
large F N , Rjj becomes independent of F N and equal to tan(r 
CUMSNR/ 2) . 

The normalized range resolution of the CW focused CLR 

measurement might be defined as the distance between two 

particular normalized percentile ranges Rj^. For example, the 
normalized 50% range resolution (dimensionless) would be 

W n (50%) £ R N (75%) - R N (25%) (28) 

and the normalized 90% range resolution would be 
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(29) 


W n (90%) ^ Rjj (95%) - Rjj(5%) . 

In general, we may use Eqs. (26) and (27) to calculate 


W N (1 - 2c) 


2 

F 1 

1 i N 

tan2$ 



L N f n J 


tan^ 

i + F 

r N - 



(30) 


where the angle 6 [rad] is 


6 - e(x - tan -1 F N ) (31) 

and 0 < e < 0.5 [dimensionless]. The normalized range resolution 
for any percentile of the signal is solved analytically. For 
W N (50%) , we insert e » 0.25, etc. For small values of F N , Eqs. 
(30) and (31) become 


W N (1 - 2e) 


2 F 


N 


tan (re) 


and for large F N , 


W N (1 - 20 


2 

tan(jre) * 


(32) 


(33) 


Note that W(50%) * W N (50%) [m] which equals 2F 2 R^ 1 and 21^ for 

the tight focus and collimated cases, respectively. 

As with Rjj, W N becomes independent of F N for large F^. 

2 

However, for small F^, W N is proportional to F N « Figure 3 plots 
W N (50%) , W N (60%) , Wjj( 70%) , W N (80%) , and W N (90%) vs normalized 

focal range F N . As predicted by Eq. (33) for large F^, Figure 3 
shows W (50%) = 2 when F N * 100. 


7. Tight Focus Case 

We now examine the special case of a tight focus, i.e., F„ 

N 
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« 1. In most cases we can set a * x. Under these assumptions 
Eq. (24) becomes 


CUMSNRfRjj) 


x + tan 


-1 


x - 


f n*n 

N ” *N 


Inverting Eq. (34) we find 


F„ tan ( xCUMSNR) 

V CUMSNR > - r” + tan ( tCUMSNR) 


(34) 


(35) 


We may also obtain Eq. (35) from Eqs. (26) and (27) . Letting 
CUMSNR « 0.5 we obtain R N (50%) - F N or R(50%) = ^1^(50%) * F, as 
expected. Half of the SNR is due to ranges from 0 to the focal 
range F. Similarly, we can calculate 


Using Eqs. 


R(25%) « F ( 1 - F n ) - 

R(75%) 8 F (1 + F n ) - 
(36) and (37), or Eq. (32), 



we also find 


(36) 

(37) 


W(50%) * R(75%) - R(25%) - 2FF N 




(38) 


This answer matches the results of others 3 * 1 ® who assumed tight 
focusing and calculated the distance between ranges where the SNR 
fell to half of its peak value. However, Eqs. (26) and (27) are 
valid under any focusing conditions and for any definition of the 
CW laser radar range resolution. It is interesting to note that 
for the tight focusing case W(50%) can, through the use of Eq. 
(43) , be written as 
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(39) 


* D min 

W(50%) = 2 J in 

which is just twice the Rayleigh range of the focused spot [see 
Eq. (2)]. Also, the product of W(50%) and the receiver solid 
angle is simply 2A, independent of focal range. 12 This is 
interesting since the SNR is approximately proportional to the 
product of the interaction length, W(50%) , the receiver solid 
angle, and ( 5 . 

8 . Beam Diameter 

Calculations of the behavior of the beam diameter are 
necessary for the measurement volume calculations in Sections 9 
and 10. The e -2 intensity diameter [m] of the beam at any range 


D(R) 


1 - 


-] 2 + < 


1/2 


(40) 


At the focal range, R = F or Rjj - F N , the beam diameter is 


D(F) 


4AF 

irD_ 


D F 
o 


D o F N 


(41) 


The e" 2 intensity diameter at the beam waist (peak SNR) position, 
R - Rp, is 


D <V = D min 


1 + (• 


. 1 — ; ) 2 ‘ 
f n * J 


1 / 2 ' 


(42) 


For a tight focus (F N « 1) 


D , a D F„ 
min o N 


(43) 


which is identical to D (F) since Rp « F. The tight focus beam 
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waist is simply the transmitted beam diameter multiplied by the 
normalized focal range. For the large F (F N » 1) case 


which together with Eq. (7) implies that the beam waist occurs at 
the transmitter and that the beam gets larger with increasing 
range . 


9 . Measurement Volume 

The e‘ 2 intensity diameter of the Gaussian beam at any 
range, D(R) , is given by Eq. (40) . The measurement volume V [m 3 ] 
may be found by integrating elements of volume dV = tD 2 (R)dR/4 
from a lower range limit RL to an upper range limit, RU: 



We define V(50%) as the measurement volume calculated with RL = 
R(25%) and RU = R(75%) , and therefore as the measurement volume 
from which 50% of the SNR arises. Plots of the 50% and 90% 
measurement volumes will be given in the next section. Inserting 
Eqs. (26) and (27) for RU and RL into Eq. (45) leads to an 
analytic expression for V, which is quite complicated. We can, 
however, make a rough estimate of V(50%) for the tight focusing 
case. We have shown that the 50% range resolution in this case 
is W(50%) * 2 FFjj [see Eq. (38)]. An approximate volume can be 
calculated by taking the volume element as a cylinder of this 
length and with a diameter D = D F„. Since at the extremes of 
the true volume element the beam diameter will be ^ 2 larger, this 
estimate will be less than a factor of 2 too small. The volume 
of this cylinder can then be written as 
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V (50%) 


(46) 


32 A 3 F 4 

From this expression we see that in the tightly focused case, the 
50% contributing volume is proportional to the optical wavelength 
cubed and to the fourth power of the focal range. It is 
inversely proportional to the fourth power of the laser radar 
beam diameter. This rapid decrease in V for smaller A and F 
values, and for larger D q values may be a very important 
consideration in cases where aerosol densities are low, such as 
at high altitudes. 

We find it interesting to note that, as far as we can 
determine, there exists no "natural" volume to which the 
calculated volumes can be normalized. The parameters R and F can 
be normalized to the Rayleigh range R^, whereas volumes cannot. 

10. Example Plots 

The equations and plots of normalized parameters provided so 
far will allow investigators to characterize their specific laser 
radar system over a wide span of parameter values. Familiarity 
and intuition can be gained from examining unnormalized parameter 
plots for a specific laser radar system. We will examine the 
case of the NASA-Marshall Space Flight Center (MSFC) monostatic 
CLR system 11 ' 13 that has recently made two Pacific ocean basin 
survey flights on the NASA DC-8 aircraft for the purpose of 
measuring the atmospheric aerosol backscatter coefficient /? as 
part of NASA's Global Backscatter Experiment (GLOBE) . The laser 
radar's parameters were approximately A * 9.13 fm, D Q = 6 cm, and 
F - 50 m; resulting in a Rayleigh range of 309.7 m. 

Figure 4 shows relative SNR vs range R for six settings of 
the focal range F [see Eq. (1)]. Close examination shows that 
the peak SNR occurs at the greatest range when F = R^ * 309.7 m, 
and is located at Rj^/2, as expected. The curves were normalized 
to have a value of 1 at R — 0. The normalization constant 
involves only D q and A , but not F . 
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Figure 5 plots CUMSNR vs range R for focal ranges F of 50, 
100, 200, 400, and 1000 m. Eq. (24) is employed with R = R^R^ 
and F = R R f N* T ^ ese curves correspond to five of the six curves 
in Figure 4. Even when focused at 50 m, at least 1% of the 
signal comes from ranges greater than 250 m. Therefore, a cloud 
or other object located at R < 250 m, with a backscatter 
coefficient 10 or more times greater than the aerosol backscatter 
coefficient, will cause errors of 10% or more. 

Figure 6 shows percentile ranges R(x) , for x = 5%, 25%, 50%, 
75%, and 95%, as a function of focal range F [see Eqs. (26) and 
(27)). The maximum plotted focal range of 200 m is still smaller 
than the Rayleigh range of 309.7 m. Therefore, Figure 6 
corresponds to the left half of Figure 2. Figure 7 gives the 
range resolution for the "center" 50% of the signal, W(50%) = 
R(75%) - R(25%) , and also the range resolutions W(60%) , W(70%) , 
W (80%) , and W(90%) vs focal range F. It corresponds to the left 
half of Figure 3. Figure 8 plots the same five range resolutions 
W vs diameter D Q when F =» 50 m. The largest values of W (poorest 
range resolution) occur for D Q ~ 2 cm. The existence of a 
diameter giving the poorest range resolution is perhaps not at 
first obvious, but it can be explained as follows. For very 
small diameters D Q , the beam will diffract quickly and hence as 
D q -♦ 0, the range contributing to SNR also goes to zero. For 
large D q , the beam can be focused tightly, and as D q -t oo the beam 
waist diameter approaches zero, and hence the contributing range 
again approaches zero. Between these two extremes lies an area 
where the range resolution W must be nonzero, and as a result, a 
maximum must exist. Figure 9 plots five percentile ranges vs 
diameter. The smallest percentile ranges appear monotonic for 
this excursion of diameter values, while the largest percentile 
ranges exhibit a peak. Therefore, the range resolution peaking 
observed in Figure 8 is due to the larger percentile ranges. 
Similar plots for F * 200 m show the poorest range resolution 
(largest W) occuring at D q « 5 cm. For very small values of D Q , 
both the percentile ranges R and the range resolutions W are 
proportional to D^. For large D Q , the percentile ranges R become 
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independent of D Q and converge to the value of F, and the widths 
W become proportional to D^ 2 . This is consistent with the small 
and large F N behavior discussed earlier, since that corresponds 
to large and small D Q behavior, respectively. 

Using Eq. (45) , Figure 10 presents V(50%) vs focal range F 
for diameter D q values of 2, 4, 6, 8, and 10 cm. For F * 50 m 
and D q * 6 cm, V(50%) « 10" 3 m 3 . Figure 11 shows the 
corresponding values of V(90%). For F * 50 m and D q = 6 cm, 
V(90%) « 7 10' 2 m 3 , or a factor of 70 larger than V(50%). For 
small values of F, the larger values of D Q have smaller 
measurement volumes due to greater focusing efficacy. The 
situation is reversed for large values of F; larger values of D q 
have larger measurement volumes, since the beams are 
approximately collimated. Figures 12 and 13 plot V(50%) and 
V(90%) vs diameter D q for values of focal range F of 10, 100, 
1000, and 10,000 m. The measurement volumes peak for some 
diameter D Q depending on F and A. Both V(50%) and V(90%) appear 
to peak at the same value of D Q . We suspect that the diameter 
yielding the maximum measurement volume is the same as that 
yielding the largest (poorest) range resolution (see Fig. 8) . 
For the special case of a tightly focused beam, one can make 
rough estimates of the ratio V(90%) /V(50%) . Such estimates show 
that the ratio is approximately 100. Numerical calculations bear 
this out, and in addition, show that the ratio stays in the very 
narrow range of 50-100 for all values of F and D Q shown in 
Figures 10-13. While it appears difficult to prove this 
statement analytically, it is nevertheless an interesting result. 
Figures 10-13 exhibit a large dynamic range of measurement 
volumes. 

11. Conclusions 

The CW laser radar SNR changes by only a factor of 2 due to 
beam diameter and focal range, and is proportional to j3A 3 . 
General analytic expressions were derived for the cumulative SNR 
vs range, for the range to a given value of cumulative SNR, for 
the range resolution of the measurement given any definition of 
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it in terms of percentile ranges, and for the measurement volume. 
The equations are valid for any parameters of the CW laser radar. 
Plots were presented in normalized parameters to cover a large 
portion of parameter space, and in unnormalized parameters for a 
specific example to gain familiarity with, and intuition from the 
curves. The behavior of the measurement volume vs focal range 
and beam diameter was plotted. The ratio of the volume 
containing 90% of the signal to the volume containing 50% stays 
in the narrow span of 50 to 100. Both the range resolution and 
the measurement volume peak for a certain diameter value. The 
measurement volumes vary over a large dynamic range. The 
question of whether small or large values of range resolution and 
measurement volume are desired depends on the specific 
application of the focused CW laser radar system. 
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14 . Figure Captions 

Figure 1. Cumulative SNR vs normalized range as a function of 
normalized focal range. 

Figure 2. Normalized percentile ranges vs normalized focal range. 

Note that R N (95%) is divided by 10 and shown as a 

dashed line. 

Figure 3. Normalized range resolutions vs normalized focal range. 

Figure 4. Relative SNR vs range as a function of focal range. 
The SNR is normalized to 1 at R * 0. 

Figure 5. Cumulative SNR vs range as a function of focal range. 

Figure 6. Percentile ranges vs focal range. 

Figure 7. Range resolution widths vs focal range. 

Figure 8. Range resolution widths vs beam diameter. 

Figure 9. Percentile ranges vs beam diameter. 
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Continuous wave single- longitudinal- mode power of 58 
mW and frequency tuning of ~1 THz has been obtained 
near 2.1 fim using a diode laser pumped thermoelectrically 
cooled Tm:Ho:YAG laser. Key words.* solid state lasers , 
single-frequency lasers , holmium lasers. 

Recent demonstrations of cw diode laser pumped thulium 
(Tm) and holmium (Ho) doped solid state lasers operating at 
wavelengths near 2 pm have shown that these lasers can 
operate efficiently at or near room temperature. 1-6 Fan et 
al. 2J first demonstrated cw room temperature operation of a 
diode laser pumped Tm:Ho:YAG laser, obtaining laser 
threshold with only 4.4 m W of absorbed pump power. Kintz 
et al. 4 demonstrated a slope efficiency of 56% in a room 
temperature cw diode laser pumped TnuYAG laser and 
Kane and Wallace 6 demonstrated a slope efficiency of 58% in 
a cw diode laser pumped Tm:Ho:YAG laser cooled to — 55*C. 
The long lifetime of the upper laser levels in TmiYAG and 
Tm:Ho:YAG lasers make efficient energy storage for pulsed 
operation a possibility. Their emission wavelengths near 
2 Mm are eyesafe, making their use in applications requiring 
atmospheric propagation attractive. In this Letter we de- 
scribe a tunable cw single- longitudinal -mode (SLM) 
Tm:Ho:YAG laser which is being developed as part of a 2.1- 
Mm coherent laser radar system. This is, to our knowledge, 
the first demonstration of SLM operation from a diode laser 
pumped TsuHo:YAG laser. 

The resonator configuration is shown in Fig. 1. The flat 
resonator mirror, Mu is a high reflector at wavelengths near 
2100 nm and transmits 88% of the diode laser pump light 
near 785 nm. Resonator mirror Ms serves as the output 
coupler transmitting ~2J>% of the laser light near 2100 nm. 
The radius of curvature of this mirror is 7.5 cm and the two 
resonator mirrors are in a near hemispheric configuration 
being separated by ~7J> cm. The CnTnuHo:YAG crystal 
used in these experiments has nominal dopant concentra- 
tions of 7.9 X 10 19 -cm~ 8 Cr 2 *, 8 J2 X XO^-cm' 3 Tm 2 *, and 5.0 
X 10 l9 -cm“ s Ho 3 *. The Cr ions are pr es e nt for fiashlamp 
pumping and do not play an active role in the energy transfer 
under diode pumping. 2 The crystal dimensions are 4 mm 
wide X 2 mm high X 3 mm long. The endfaces of the crystal 
are polished fiat and AR coated with the power reflectivity 
being <0.2% at wavelengths near 2100 nm and <5.0% at the 
785-nm pump wavelength. The 3- mm long laser crystal 
transmits 28% of the incident 785- nm diode laser pump light 
The crystal is mounted on a three-stage thermoelectric (TE) 
cooler allowing its temperature to be varied from +30 to 
-40*C. All the data reported below is for operation at 
— 40*C. The TE cooler and laser crystal are enclosed in a 
nitrogen-purged chamber and one end of the laser crystal is 
~0.5 mm from Ml. 

The uncoated solid fused silica etalons, E\ and have 
thicknesses of 1.0 and 0.09 mm, respectively. The thin eta- 
Ion, Eu is used to control the 1 aeer wavelength and the thick 
etalon, E\ f is used to obtain SLM operation. The fused silica 
Brewster plate controls the polarization of the laser. The 
laser is longitudinally pumped with the output from a ten- 


Brewster 

□i a f e Cr : Tm : Ho : YAG 




Fig. 2. Laser output power near 2100 am v» diode User pump power 
near 785 nm incident on rear cavity mirror A/,. 


stripe 500-mW cw phase coupled diode laser array (Spectra 
Diode Labs SDL-2432-P1). The diode laser wavelength is 
tuned to the TnuHoiYAG absorption maximum near 785 nm 
and controlled there using another TE cooler which is inte- 
gral to the diode laser package. The optical train that cou- 
ples the diode laser output into the laser crystal delivers 
~85% of the diode laser output to the resonator end minor, 
Mt, and focuses it to e spot size of ~68 X 37 pm in the 
CnTnrHo: YAG crystal. 

With the two intracavity etalons removed, the output from 
the laser is made up of several longitudinal modes at both 
2090 and 2096 nm. The output power vs pump power inci- 
dent on M\ for this TEMoo multiple-longitudinal -mode 
(MLM) operation at — 40 # C is shown in Fig. 2 (squares). 
The oscillation threshold occurs at 100 mW and the slope 
efficiency for points well above threshold is 33%. The non- 
linear behavior near threshold is characteristic of quasi - 
three-level laser operation. 6 * 7 Approximately 65% of the 
pump power incident on Mi is absorbed by the laser crystal. 
Therefore, a plot of output power vs absorbed pump power 
would yield a slope efficiency of 51%. This compares favor- 
ably with the theoretical maximum slope efficiency for 
Tm:Ho:YAG lasers of 75% (Le., g, - 2*^ where n* “ A^/Xl - 
785/2090 “ 37.6% end the factor of 2 it for perfect twofor- 
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Fig. 3. Output power vs wavelength obtained by tilt tuning the 
0.09-mm thick intracavity e talon Ez and with the 1.0-mm thick 
intracavity etalon E[ removed from the cavity. 


one cross relaxation of the Tm ions 3 ). Due to the Brewster 
plate and a slight tilt of the laser crystal in the near hemi- 
spheric resonator, the TEMoo laser output is slightly astig- 
matic. The degree of astigmatism is a strong function of the 
exact mirror separation. When the SLM data in Figs. 1-4 
were taken, the ellipticity was —1.16 X 1.0, and the TEMoo 
spot size on M\ (1/e 2 intensity radius) calculated from the far 
field divergence was —81 X 70 pm. 

If the thin (0.09-mm thick) etalon is inserted into the 
cavity, the laser typically operates on three to four non&dja- 
cent longitudinal modes spread out over 30 GHz. The laser 
wavelength depends on the incident angle between the axis 
of the laser mode and the normal to the etalon surface. 
Figure 3 is a plot of output power vs wavelength obtained by 
angle tuning the thin etalon. The pump power incident on 
Mi during the collection of data for Fig. 3 was 430 mW. The 
initial tilt angle in Fig. 3 is -5° at X * 2100 nm and increases 
as X decreases. The lower output power at 2090 nm com- 
pared with 2097 nm is due to increased walk-off loss due to 
the higher tilt angle of the thin etalon. Operation with 78 
m W of output power near 2090 nm can also be obtained with 
a tilt angle of —2°. 

With all the other elements shown in Fig. 1 in place, 
insertion of thick etalon E\ into the cavity at near normal 
incidence results in SLM operation as shown in Fig. 4. Fig- 
ure 4(A) shows the laser spectrum as measured with a scan- 
ning Fabry-Perot interferometer. The two peaks are sepa- 
rated by the 15-GHz free spectral range (FSR) of the 
interferometer and the linewidth is due to instrument resolu- 
tion. The longitudinal modes of the laser are separated by 
—2 GHz and, if other longitudinal modes were present, they 
would easily be resolved by the Fabry-Perot Figure 4(B) 
shows the spectrum as measured by a scanning confocal 
interferometer having a FSR of 300 MHz. Again the 1.7- 
MHz FWHM linewidth is due to instrument resolution. 
The laser was operating at 2090.3 nm, as measured with a 
0.25-m monochromator, and the output power was 55 mW 
when the spectra of Fig. 4 were recorded. The output power 
vs pump power incident on M\ for SLM operation at 2090 nm 
and at — 40°C is shown in Fig. 2 (triangles). Oscillation 
threshold occurs at 162 mW and the slope efficiency for 
points well above threshold is 27%. By angle tuning thin 
etalon SLM operation was demonstrated at several 
points in the tuning range between 2088 and 2099 nm with a 



TIME (2 nu*c/DIV) [FREQUENCY (3.14 GHz/DlV)] 



TIME (0.2 mwc/DIV) [FREQUENCY (1.09 MHx/DIV)] 


Fig. 4. Spectrum of the single- longitudinal -mode Tm:Ho:YAG la- 
ser. (A) The spectrum as measured using a scanning Fabry-Perot 
interferometer. The peaks are separated by the 15-GHz free Bpec- 
tral range of the interferometer and the linewidth is due to instru- 
ment resolution. (B) The spectrum as measured with a scanning 
confocal interferometer whose free spectral range is 300 MHz. The 
linewidth is again due to instrument resolution. 


response vs wavelength very similar in shape to that shown in 
Fig. 3. 

If the 2.5% output coupler (Mj in Fig. 1) is replaced with a 
0.5% output coupler, the laser can also be made to operate 
SLM and tune at wavelengths near 2121 nm. This depen- 
dence of operating wavelength on output coupling in 
Tm:Ho:YAG lasers has been described by Fan et al . 3 

In summary, we have achieved broadly tunable SLM oper- 
ation from a TE-cooled diode laser pumped Tm:Ho:YAG 
laser operating at -40 d C. No special precautions have been 
taken for long term stability with the current hardware. The 
laser mirrors are mounted on an aluminum baseplate and the 
TE cooler for the laser crystal is not actively stabilized. 
Even so, the laser runs for many tens of minutes before two- 
mode operation or longitudinal mode hops occur. In the 
near future, we will replace the aluminum baseplate with a 
super Invar baseplate and actively stabilize the TE cooler to 
achieve long term stability. 

This research was part of a Small Business Innovation 
Research program monitored by the Defense Meteorological 
System’s Advanced Technology Office at the U.S. Air Force 
Space Division, Los Angeles Air Force Base, California. 
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ABSTRACT 

Sirgle- frequency Tm.Ho: YAG and TnsYAG lasers 
operatlre near 2 fm are attractive sources 

for several applications including eyesafe 
laser radar (Ildar) and puplig of AgGaSes 
parametric oscillators for efficient 
generation of longer wavelengths. As part 
of a program to develop a coherent lidar 
system using TYb.Ho: YAG lasers, we have 
developed a diode laser-pimped, tunable, CW, 
single- longitudinal-mode (SIN) Cr: Tv Ho: YAG 
laser and a flashlamp-pmped, single- 

transverse- mode, Q-switched Cr,Tm,Ho: YAG 

laser. The CH laser has been used to 
injection seed the flashlanp-pmped laser, 

resulting ir SIN Q-switched output. 

Operational characteristics of the CH and 
Q-switched lasers and injection-seeding 
results are reported here. 


SLAVS CBCHJjATOR 


The f las hlanp- pimped Q-switched slave 
oscillator utilizes a polarization-coupled , 
large- volume-TEMoe cavity configuration. 
(See Slave Ctecillator in Figure 4. ) The end 
mirrors are flat high reflectors at 2. 09 (m 


airi the cavity is dynamically stabilized 
using the thermal lens of the pimped laser 
rod. The YAG laser crystal has nominal 
dopant concentrations of 7.9 x 10 1# cm" 
Cr* , 8.2 x 10” cm"* Tm*, and 5.0 x 10 1# 

cm - * Ho**. In order to produce the long 
duration Q-switched pulses (1. e. , the narrow 
Fourier- transform- 11ml ted bandwidth) re- 

quired for the coherent lidar application, 
the cavity mirrors are separated by ~ 1.2 m. 
The output energy vs input energy for 
operation at ~ 20° C is shown in Figure 1 
for both normal mode and single- pulse 


Q-switched operation at a HU? of ~ 3 Hz. 
Note that ~ 5 OX of the normal node energy is 
extracted in single Q-switched pulses. This 
is expected for the dopant concentration 
utilized since the active Ho ions only store 
~ 5 OX of the energy near room tenperature, 
with the Tm ions storing the other 5 OX. The 
transfer time from Tm to Ho is estimated to 
be > 5 fB, which is nuch longer than the 
Q-switched pulse le%ths. This slow trans- 
fer can result in multiple pulses when the 
laser is Q-switched; 1 to avoid this we gate 
the voltage applied to the Ebchels cell, 
bringing the laser cavity beck to a low-Q 
state before the second pulse evolves. 

The low slope efficiency shown in Figure 
1 is due to the combination of an 
Inefficient flashlaip pimp chamber and a 
very lossy Ebchels cell (bad AR coating). 
The 5. IX slope efficiency demonstrated by 
Quarles et al. , 2 in a multiple transverse 
and normal mode Cr, Tm.Ho: YAG laser, leads us 
to expect an Increase of a factor of four or 
more in slope efficiency when the cavity 
losses are minimized and a more efficient 
pimp chamber is utilized. The pulse 
duration for 30 mJ and 50 mJ Q-switched 
pulses Is — 250 ns and w 150 ns FWHM, 
respectively. Although no damage has been 
experienced, the Q-switched pulse energy has 
been Hmtterf to M 50 nJ to date for fear of 
damaging intracavity optics. 

At the 3 Ik FRF used throughout this 
study, the thermal lens created by the rod 
is correct for large-volume TEMoo operation 
and excellent mode quality Is obtained. The 
intensity profile of the slave oscillator 
output, measured by scanning a 100- jm 
diameter lor^- wavelength InGaAs photo- 
detector across the beam, is shown in Figure 
2. Also shown is a Gaussian fit to the 
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INTENSITY PROFILE 



Figure 1. Output pulse energy for TEMoo Figure 2. Intensity profile of the 

operation vs energy supplied to the flash- TEMoo output froa the slave ocsillator. 

lamp for both normal-mode and single- pulse 
Q-switcbed ( c ir c les) operation. 



Figure 3. Schematic of the optical configuration used in the 
s ingle- frequency raster oscillator. 
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data. Operation of the slave laser in any 
one of the four wavelength regions near 
2090, 2096 , 2121, and 2128 im can be 
obtained by tilt- tuning the intracavity 
thin-film polarizer. (The transmission of 
the polarizer for p- polarization light 
varies with wavelength and the spectral 
transmission peak can be tuned by small 
changes in the tilt angle. ) 

MASTER OSCILLATOR 

The diode lasez^ pumped, CW master oscillator 
3 town in Figure 3 utilizes a linear cavity 
with an intracavity Brewster plate for 
polarization control, and two solid 
f used-silica etalons (thicknesses of 0.5 and 
0. 09 mm) for frequency control. The 0. 09 ran 
thick etalon, Ea, is unooated and the 0. 5 ran 
thick etalon, Et , is coated for 18X 
reflectivity at 2. 09 pm. The intracavity 
lens, Li , (focal length = 75 mm) is utilized 
to provide stability and to minimize the 
waltoff loss introduced as the 0. 5 ran thick 
etalon is tilted. Both mirrors , Hi and Ha, 
axe piano /piano with Hi highly reflecting 
and Ha ~ 97. 5% reflecting at 2. 09 pm. The 
entire optical train is mounted in a 
t emp e r a t ure con t rol lad super invar block. 
Other than higher frequency stability, the 
other operational characteristics of this 
laser are virtually identical to those 
described in reference 3, som only provide 
a summary here. With the laser crystal 
thermo-electrical ly cooled to ** 230 K, to 
have obtained M 50 nH of SU4 output power 
near 2090 m and 2096 nm using 430 nH of 
pump power from the diode laser. A scanning 
oonfocal interferometer hem shown that the 
linewidth and short-term frequency jitter 
from this laser are less than 1 tffe. By 
tilt- tuning the intracavity etalons, tuning 
over the wavelength range between 2067 and 
2099 nm has been achieved. Tuning near 2121 
nm has also been demonstrated, using the 
laser described in reference 3, fay changing 
the reflectivity of the output coupler. 

INJECTION SEEDING 

The apparatus used for the injection- 
seeding studies is illustrated in Figure 4. 
The output of the master oscillator was mode 
matched and injected into the slave 
oscillator using the "back" of the 
Intracavity thin-film polarizer. Isolation 
of £ 40 dB is provided by a yittrium iron 
garnet Faraday Isolator. Neutral density 
filters can be placed in between the master 
and slave oscillators to reduce the injected 
intensity. The output of either laser can 
be analyzed using seaming interferometers 
and a quarter-meter monochromator. A fast 


long-wavelength InGaAs photodetector is used 
to observe the tenporal profile of the 
pulsed output. The longitudinal mode 
frequency of the slave laser is controlled 
by translation of a slave cavity end mirror 
which is mounted on a piezoelectric lusher. 
Throughout this study the required match 
between the master oscillator frequency and 
the slave cavity longitudinal mode frequency 
was maintained by manually controlling the 
voltage applied to the end mirror P27T. The 
"man in the loop" adjusted the PZT voltage 
so as to minimize the observed pulse buildup 
time. 4 A second quarterrove plate can be 
added to the slave oscillator (shown as the 
optional A/4 plate in Figure 4) to minimize 
spatial hole burning in the laser rod. 5 For 
the results presented here, the optional 
waveplate was not utilized. 

Figure 5 shows an example of the temporal 
evolution of the Q-switched slave laser 
pulse for both seeded and unseeded operation 
with the wavelength near 2090 nm. The 
right-most 30 mJ unseeded pulse exhibits 
longitudinal mode beating at ~ 120 MHz 
(shown as M 20 MBs due to aliasing by the 
100 magasanple-per^-seoond digital oscil- 
loscope used). Note the earlier buildup 
time and smooth profile for the left-most 
seeded pulse. The lack of mode beating on 
the seeded pulse suggests SUi output. In 
this exanple, the seeded laser pilse 
contains « 7 OX of the energy that the 
unseeded laser pulse contains. This 
reduction is probably due to a combination 
of spatial hole burning in the laser crystal 
and the lower gain experienced due to the 
detunlre of the seed laser frequency from 
line ooiter (see next paragraph). 

With the experimental conditions the same 
as those used to produce Figure 5, a 
scanning Fabry Perot (FP) etalon whose free 
spectral rar«e (FSR) is 75 (Ms was used to 
measure the spectrum of the slave and master 
oscillator outputs (see Figure 4). Figure 6 
shows the spectrum measured during a 50 
second scan with the master oscillator 
output blocked from the slave oscillator. 
The "grass" to the left is due to the 
unseeded Q-switched slave laser output and 
the smooth peak to the right is due to the 
master oscillator output. The width of the 
smooth peaks represents the resolution of 
the interferometer and the repetitive nature 
of the display results from scanning by more 
than 1 FSR. The right-most feature on the 
display is due to resetting the FP mirror 
spacing to begin another scan. Figure 6 
shows that for this example the bandwidth of 
the unseeded Q-switched output is ^ 10 GHz 
and the frequency of the CW master 
oscillator is offset by ~ 22 CSfe from the 
centroid of the slave laser output 
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Figure 4. Experimental apparatus used to study injection seeding. 



TIME [200 nS/DIV] 


Figure 5. Temporal evolution of the unseeded (right) 
and seeded (left) Q-switehed pulses. 
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frequency. If the output of the CW master 
oscillator is allowed to enter the slave 
oscillator and if the seed radiation is 
resonant or nearly resonant with a 
longitudinal mode frequency of the slave 
oscillator, the spectrum of the slave laser 
output collapses and is pulled to the 
lo ngitudinal mode which is closest to that 
of the seed radiation. This is illustrated 
in Figure 7 where the experimental 
conditions were identical to those described 
in the discussion of Figure 6, except that 
tte master oscillator was allowed to seed 
the slave oscillator. The peak voltage seen 
in Figure 7 represents the saturation 
voltage of the detector circuitry. 

Higher resolution spectra of the pulsed 
output are shown in Figures 8 and 9 for both 
unseeded and seeded operation. For this 
measurement a scanning oonfocal interfero- 
meter whose FSR is 300 Hk was used. The 
upper trace in both figures shows the 
50- second voltage ramp used to drive the 
scanning interferometer. Note that when 
seeded the spectrta collapses to < 6 Mfe 
FWW of which a large part is due to the 
resolution of the interferometer. The 
transmitted intensity scale in Figure 9 is 
ten times larger than that of Figure 8. The 
spectrum of Figure 8 is uniformly filled due 
to alleging of the 10-20 Ofe wide unseeded 
spectra into the 300 Mfe FSR of the 
interferometer. 

0y tuning the wavelength of the master 
oscillator, we have achieved SU4 operation 
over a *» 250 GHz (3. 7 nm) wide band centered 
at 2090 nm and ~ 100 CHs (1.5 m) wide bands 
centered at 2096.2 nm and 2097.7 nm. At the 
center of each of these bands, the energy of 
the seeded SW pulses is typically 80-90X of 
that of the u n soo d od pulses. At the edges 
of each wavelength band, the seeded energy 
Is n 5 OX of the unseeded energy. Seeding 
near 2121 nm and 2128 nm has not been 
attempted at this time. 

A measurement of the maxima detuning 
between the frequency of the nmster 
oscillator and that of the nearest slave 
cavity longitudinal mode, which results in 
SLH operation, was wade as a function of the 
injected power. This maxim a detuning is 
ooamonly called the injection-seeding (or 
lockii«) range. 8 A large injection-seeding 
range is useful since it eases the 
requirements an the servo system locking the 
frequency of the master oscillator to a 
slave cavity longitudinal mode. The 
injection-seeding range vs injected power is 
shorn in Figure 10 far seeding at 2097.8 in 
when the unseeded pulse energy was ~ 25 mJ 
and the seeded pulse energy was M 20 ad. 
Note that with only 1 nW of injected power 
the injection-seeding range is ~ 10 Mfe. He 


have also obtained SUM Q-switched operation 
with only 670 nW of injected power; this 
point is not shown because the injection- 
seeding range was so small (< 2 Mfe) that a 
reliable measurement could not be easily 
made. 

SUIWARY 

He have demonstrated tunable injection- 
seeded SLH operation of a room tenperature 
Q-switched Cr,Tta,Ho: YAG laser. To 

aoconplish this task a highly-stable tunable 
s ingle- frequency CW master oscillator was 
constructed and a tunable TEMoo Q-switched 
flashlanp-pumped oscillator which produces 
over 50 mJ/pulse was constructed. In the 
near future we will construct electronics to 
automatically servo the slave cavity 
longitudinal mode frequency resulting in 
unattended SUf Q-switched operation of the 
slave laser. 
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Figure- 6. Qnseeded spectra as observed Figure 7. Seeded spectra as observed using 
usiiH a seaming interferometer whose free a scanning interferometer whose free 
spectral range is 75 GHz. spectral range is 75 Ufa. 




Figure 8. Qnseeded spectra as aeesured Figure 9. Seeded spectra as measured using 
teli*! a s eamin g interferometer whose free a scanning interferometer whose free 
spectral range is 300 Ufa. spectral range is 300 Ufa. The transmitted 

intensity scale is ten times that of 
Figure 8. 
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Figure 10. Injection-seeding range vs injected power far operation near 2097.8 not 
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APPENDIX G 


0 CALIBRATION CONSIDERATIONS FOR THE 
AFGL PULSED C0 2 COHERENT LASER RADAR 

Michael J. Kavaya 
Coherent Technologies, Inc. 

An AFGL ground-based pulsed C0 2 coherent laser radar (lidar) 
has been used and will be used in the future to obtain calibrated 
profiles of the atmospheric aerosol backscatter coefficient 0 
(Alejandro et al., 1990). (See Section 11 for references.) This 
appendix attempts to list the legion of calibration choices and 
trade-offs faced by the lidar designer and operator. The typical 

trade-offs for improved calibration are: increased cost, 

increased complexity of lidar design, increased difficulty, 
and/or increased measurement/calibration time. The overall goals 
are: 

1. Reduce random noise 

2 . Model measurement physics and reduce systematic errors 

3. Obtain accurate system zero 

4. Obtain accurate system gain 

1. Telescope Overlap Function 

The telescope overlap function (TOF) , 0(R) , also referred to 
as crossover function or geometrical form factor, is a desired 
term in the coherent lidar equation, for received signal power or 
SNR, that is used to process received data and produce a y0(R) or 
0(z) profile. Ideally, 0(R) represents all range dependent 
effects on the signal that have not been included in the lidar 
equation and data processing steps. In this perspective, the 
name telescope overlap function is a misnomer. The lidar 

equation employed, P r (R) or SNR(R) , may originate from analytic 
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theory or from computer numerical analysis, or may combine one of 
these techniques with a heuristic modification of the equation. 
A past analysis of 0(R) for a specific C0 2 coherent lidar system 
(Kavaya and Menzies, 1985) revealed that incorrect use or neglect 
of 0(R) could result in values of P in error by factors up to 30. 
Since that study was performed, newer lidar equations have become 
available (Frehlich and Kavaya, 1990) that include more of the 
range-dependent physics of the coherent lidar measurement. 
Therefore, the TOF is an evolving function, and is specific to 
each lidar system. It may or may not need to include the effects 
of misalignment of the transmitted (XMTR) beam and the back- 
propagated local oscillator (BPLO) beam/ a bistatic lidar with 
separated XMTR and BPLO; beam truncation by optical elements; 
obscuration by the telescope (if any); the spatial profile of the 
XMTR or BPLO beams, if not Gaussian; the detector size, shape, 
and spatial profile of quantum efficiency; mis-positioning of the 
detector (related to misalignment) ; and pulsed profile temporal 
shape effects. For a bistatic lidar measuring at short ranges, 
0(R) will be significant due to the XMTR/ BPLO beam separation. 

The best determination of 0(R) is to fire the lidar into an 
unchanging atmosphere with well known total extinction oc(R) and 
aerosol backscatter $(R) profiles. With enough shots averaged to 
reduce the effects of speckle; the SNR(R) , /?(R) , and a(R) 
profiles, and the chosen lidar equation could be used to 
calculate 9 ( R) . This may be impossible for the present. A 
relative (vs absolute) 0(R) curve may be determined if the lidar 
is fired into a homogeneous atmosphere where /?(R) = P and a(R) = 
a, both constants vs range. The most likely scenario is to model 
0(R) to include as much physics as possible that is not included 
in the lidar equation (Ancellet et al., 1986). This absolute 
model of 0(R) may then be checked by varying the settings of the 
lidar and comparing experiment to theory, by firing into an 
assumed uniform atmosphere and looking for a uniform P(R) result, 
or by varying the elevation angle of the lidar and comparing the 
resultant P values at equal heights. 
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Calibration Target 

The method of choice for calibrating pulsed coherent lidar 
measurements of 0 is to use a hard calibration target (Kavaya et 
al., 1983; Kavaya, 1985; Kavaya and Menzies, 1985; Kavaya, 1987). 
The use of the total backscatter to Rayleigh backscatter ratio 
calibration technique, which is used in the visible wavelength 
region (often in concert with a hard target anyway) , is not 
feasible at CO 2 wavelengths since molecular backscatter is 
negligible. The calibration technique that assumes a functional 
relationship between backscatter and extinction (a field of 
research in itself) does not apply since CO 2 backscatter is 
dominated by aerosol particles while extinction is dominated by 
molecules. The rule of thumb that says to calibrate an 
instrument with a known target as close in character as possible 
to the eventual target (e.g., in this case, aerosol particles, 
oil drops, homogeneous spheres, etc.) is not practical since it 
is difficult to make a well characterized target that mimics 
aerosol particles, and since this calibration target would have 
to fill the expected measurement volume, which is more than 150 m 
long for a l-/xs long pulse. The technique of using a 
comprehensive lidar equation and then characterizing the 
reflectance or transmittance or gain of every component is very 
hard and prone to errors as components change with time. 

The ideal calibration target would be easy to fabricate, 
inexpensive, lightweight, durable, and be as close to diffuse or 
Lambertian in its scattering behavior as possible. A large 
target area is desired so that it may be turned with respect to 
the lidar beam (to avoid specular reflectance) while still 
projecting a sufficient area to encompass the entire beam. The 
target should be turned about a vertical axis to avoid 
backscatter from the ground. While the field target should be 
large, a small witness sample with exactly the same reflectance 
behavior is desired to allow laboratory calibration of the 
desired reflectance. The laboratory calibration of the witness 
sample should employ conditions that match, as closely as 
possible, the use of the large field target in calibrating the 
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lidar system. Specifically, continuity should be sought for: 

1. the target itself, its reflectance characteristics, 

2. the laser wavelength and linewidth, 

3. the incident and reflected (detected) polarization 
states, 

4. the polar angle of the laser beam with respect to the 
normal to the target (45° is a nice compromise between 
avoiding specular reflection and losing projected area, 
i.e., projected area = 71% of actual area), 

5. the incident and detected solid angles (as close as 
possible) , 

6. the size relationship of the target area, the 
illumination area, and the detector's field of view (it 
is recommended that the target's projected area be the 
largest, and be sufficiently big that the beam is 
always contained in the projected area despite laser 
pointing jitter and beam wander due to atmospheric 
refractive turbulence) , 

7. complete averaging of speckle. 

The target should exhibit no opposition effect (cube corner 
behavior) since the necessarily different solid angles employed 
in the laboratory calibration of the witness sample and the field 
calibration of the lidar will lead to errors in the target 
reflectance parameter p* , and hence errors in /?. 

It is desirable to ensure that the target is moving randomly 
over a distance of many wavelengths of light in order to average 
speckle effects. Often the wind or vibration from traffic will 
guarantee this. The probability distribution function (PDF) of 
the received IF power from the target shots should be checked to 
confirm negative exponential behavior. 

The optical train used for atmospheric and hard target 
measurements should be identical. Any differences should be 
compensated for. 

For CO 2 laser work, the researcher should employ the Lidar 
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Target Calibration Facility at NASA's Jet Propulsion Laboratory 
(Kavaya, 1985) . Presently, there is no charge for calibration of 
witness samples. 

3 . Atmospheric Extinction 

Proper derivation of the /?(R) profile requires knowledge of 
the total extinction profile a(R) that was encountered by the 
lidar beam. Components of extinction include molecular 
absorption, molecular scattering, aerosol absorption, aerosol 
scattering, and the effects of clouds. Ideally, a(R) would be 
measured independently for each lidar shot and used in data 
reduction. Practically, however, it is difficult to obtain oc(R) 
even as an average over the measurement time. Models of the 
atmosphere may be employed to calculate oc(R) for certain assumed 
conditions. This is probably the only practical choice at 
present. However, the atmosphere models do not account for a 
ground-based lidar shooting through a boundary layer. It is 
highly recommended that a separate measurement run of the lidar 
(aiming slightly above horizontal) be employed to determine the 
total extinction coefficient (assumed constant) and height of the 
boundary layer (Kavaya and Menzies, 1985) . 

4. Lidar System Zero 

If one could order an atmosphere with 0 = 0, then the range 
dependent system zero profile or minimum detectable 0 profile 
could be experimentally determined. Since this is impossible, 
one has to choose a method of determining this curve. This 
information is needed to subtract any DC offsets from the data 
before processing, and to display a minimum detectable 0 curve 
along with the calculated 0 curve. 

Choices available include blocking the transmitter 
telescope, using pretrigger recorded data (before the laser 
fires) , or using data at ranges so large that there can be no 
signal. Blocking the telescope may create scattered light and 
electronics saturation that is not present during actual 
measurements. The pretrigger data may not include noise caused 
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by the firing of the laser. The large range data may miss 
ringing of the electronics caused by firing the laser. The 
researcher should examine signals from all three techniques to 
determine the best approach. 

5. Effects of Transmitted Pulse Profile 

The effects of the actual temporal profile of the 
transmitted laser pulse can affect the accuracy of the calculated 
/?( R) profile (Kavaya and Menzies, 1985) . The most accurate 
approach is to employ a lidar equation that has not assumed 
anything about the pulse shape or duration. However, if the 
pulse shape and duration are repeatable, a range-dependent 
correction factor may be calculated and applied during data 
reduction. Unfortunately, this correction profile depends on the 
total extinction profile, oc(R). 

6. Averaging a Nonlinear Function of Received Power 

Because of speckle fluctuations, many lidar shots are 
averaged to obtain an estimate of the mean profile of received IF 
power. However, care must be taken since often the voltage into 
the A/D converter is a nonlinear function of IF power. Examples 
include outputs voltages that are the square root of power or the 
logarithm of power. Averaging these voltages over many pulses 
and performing the inverse operation on the average, e.g., V 2 or 

10 v , will cause errors in /?(R) (Kavaya and Menzies, 1985) . 
Ideally, it would be arranged to average a voltage linearly 
proportional to IF power. If more dynamic range is needed, then 
parallel A/D converters, or a switch to higher gain during the 
return signal, or a ramping up of gain should be considered. If 
this is not possible, then correction profiles for /?(R) can be 
calculated from knowledge of the nonlinear relationship and the 
PDF of the detected voltage. 

7. Pulse Energy Normalization 

Ideally, the received signal from every shot would be 
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normalized by the pulse energy of that shot. Summing the pulse 
energies independently and dividing the summed signals by the 
energy sum may cause systematic error. If the laser is known to 
have sufficiently small fluctuations in transmitted energy, then 
the normalization step may be skipped, and appropriate error bars 
added to the /?(R) profile. 

8. Gain Reduction for Hard Target Data 

The reflectance p* [sr _1 ] of hard targets is typically so 
large compared to atmospheric (actually 0ct^/2 sr" 1 ) that the 
lidar receiver gain must be decreased. Changing the telescope 
aperture in any way is not advisable since this changes the 
physics of the heterodyne detection, and hence changes the 
correct lidar equation. A convenient technique is to reduce the 
receiver gain by a known ratio (i.e., lower the gain by a known 
number of dB) . Preferably, this gain reduction is accomplished 
with optical neutral density filters or low reflectance 
beamsplitters rather than with lower electronic gain, so that 
target calibration is performed under conditions similar to high 
/? measurements. (The larger output of the beamsplitter could be 
used for a direct detection measurement for estimating p*.) This 
technique is fine provided the recorded voltage is proportional 
to the IF signal power or its square root. Then the ratio of the 
two gains employed is all that is required. However, the use of 
a logarithmic amplifier in conjunction with this gain ratio 
technique can make proper calibration very hard (Kavaya and 
Menzies, 1985) . The absolute values of the two receiver gains 
are then needed, which are much harder to obtain than their 
ratio. 

9. Receiver Gain vs Signal Size 

The lidar receiver should be mapped out to determine if its 
gain varies with signal size, especially for low signals. If the 
receiver gain is not constant, then the predetermined map should 
be used during data reduction. 
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10. Receiver Gain vs Signal Frequency 

The return signal from a hard target duplicates the temporal 
profile of the transmitted pulse. The atmospheric return signal 
varies much slower in a mean value sense, but does have speckle 
fluctuations imposed on it. Care should be exercised to ensure 
the effect of the spectrum of the received signal on the receiver 
gain is understood. 

11. Receiver Saturation 

Saturation of any of the receiver components will cause 
errors in /?(R) . The absence of saturation in each component 
should be checked. 

12. Detector Operating Point 

The DC and AC load lines imposed on the photodetector by the 
bias circuit and preamplifier should be chosen to provide linear 
operation over the required dynamic range of the signal (Post and 
Cupp, 1990) . 

13 . Shot-Noise-Limited Operation 

The term shot-noise-limited is ambiguous. Most researchers 
simply mean that shot noise is some factor larger than all other 
sources of noise together. The adjective has also been misused 
to mean that shot noise equals the sum of other noise sources. 
It should be arranged that shot noise is a least 10 times larger 
than the other noise. Otherwise, a correction factor should be 
inserted into the lidar equation (Post and Cupp, 1990) . 

14. Detector Quantum Efficiency 

Many HgCdTe detectors are delivered with only the DC quantum 
efficiency specified. The actual AC quantum efficiency at the 
signal IF frequency should be employed in the lidar equation. 

15. Misalignment 

Both the lidar equation and the TOF model should include 
misalignment effects (but include each effect only once) . Data 
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should then be taken vs parameters of the lidar system such as 
XMTR and BPLO radii and focal ranges, and detector (x,y,z) 
position to confirm correct alignment. Maximizing the signal 
from a hard target at some close range does not provide good 
alignment. 

16. A-Scope Capability 

The researcher should be able to observe a pulse averaged 
A-scope display to confirm lidar operation 

17. Polarization of Light 

The fact that 0 is really a 4 x 4 Mueller matrix and that 
each optical component affects the polarization of the laser 
light should be accounted for (Kavaya, 1987; Anderson, 1989) . 

18 . Neglected Processes 

The usual lidar equation neglects multiple scattering, 
Rayleigh scattering, fluorescence, resonant fluorescence, Raman 
scattering, and any effect of the beam on the atmosphere itself. 

19 . Refractive Turbulence 

Refractive turbulence effects may affect even C0 2 lidar 
systems when a long horizontal path is used for calibration. The 
lidar equation should include these effects (Frehlich and Kavaya, 
1990) . For collimated beams, the importance of refractive 
turbulence can be determined by measuring the scintillation of 
the on-axis transmitted irradiance at the calibration target. 

20. Secondary Calibration Technique 

A secondary calibration technique, that is much easier to do 
than the primary hard target technique, is desirable in order to 
provide frequent checks on lidar gain stability. 

21. Data Quality Monitor 

A real-time data quality monitor would alert the lidar 
operator of problems, and could save considerable amounts of 


120 


time. 


22. Signal Processing Hardware 

Digital signal processing may provide more accuracy and 
flexibility than analog processing (e.g., a saw processor). 

23. SNR or Power Calibration 

Calibration using the aerosol and target SNR values removes 
any gain drifts in the system provided that shot noise is 
proportional to the system gain. Calibration using IF power 
values will be more accurate if the system gain is constant since 
estimation of only S - N is needed rather than (S - N)/N. 

24. Choice of Signal Spectrum Location 

The typical heterodyne coherent lidar system has a frequency 

offset f between the pulsed and LO lasers. The nominal IF 
o 

signal frequency from a nonmoving target is therefore f Q . In 
addition, the receiver bandwidth is designed to accept IF signals 
over the span of f Q ± ^max^' where V max is the highest allowed 
unidirectional radial wind velocity, which causes a Doppler shift 
in the received optical frequency. There are three reasons to 
choose a large value of f Q . First, electronic components are 

less expensive for smaller fractional bandwidths. In this case, 
the fractional bandwidth would be (4 v max )/(^ f 0 ^* Second » ^ 
easier to achieve a flat noise floor with smaller fractional 
bandwidths. A flat noise floor is desirable for many signal 
processing steps, such as noise subtraction. Third, the receiver 
will experience direct detection current as the laser pulse 
scatters during its transmission out of the lidar, and also when 
large returns from hard targets are received. The spectrum of 
the direct detection current will be approximately the spectrum 
of the transmitted pulse profile, i.e., DC to some upper 
frequency f u » It is important that the direct detection signals 
do not saturate the receiver electronics, and therefore it is 
recommended that f « f Q - 2 v max /^* 
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